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Abstract 


This' report outlines the preliminary design and characteristics of a 
hypersonic aircraft for a flight at Mach 10 using only scramjets for two 
minutes at 100,000 feet. There are many design problems that have to be 
addressed for such a highspeed flight. These include aerodynamic, 
thermal, logistical and structural problems. 

This report contains ideas to deal with these problems that have been 
examined by our research team, the gold team. Aerodynamic calculations, 
logistical solutions are presented along with thermal and structural 
designs. Many ideas for hypersonic aircraft are based on theory and have 
limned experimental foundation. 
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NOMENCLATURE 


AR 

b 

C D 

C Do 

C Df 

Cf 

C Dw 

C Dle 

C Db 

C Dp 

^Dp 
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C DA 

C DA(NC) 

Cl. 

V 

C mac 


Aspect Ratio 

Span of Airplane ft 

A/C Drag coefficient 

Coefficient of Drag at Zero lift 

Skin Friction Drag Coefficient 
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Leading Edge Bluntness term 

Base Pressure Drag Coefficient 

Subsonic Pressure Drag Coefficient 

Supersonic Wave DragCoeff icient 

Nose Wave Drag Coefficient 
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Interference Drag Coefficient 

Lift Coefficient 

Lift Curve Slope 

Moment Coefficient About 
Aerodynamic Center 



C mcg 

Moment Coefficient About C.G 


Cn v 

Slope Of C m vs. *-Curve 

/deg 

c r 

Root Chord 

ft 

C T 

Tip Chord 

ft 

d 

Diameter of A/C 

ft 

dT 

Time of Impulse 

sec 

e 

Wing Efficiency Factor 


e' 

Wing Planform Efficiency Factor 


°F 

Temperature 

degrees 

f N 

Nose Fineness Ratio 


f A 

After-body Fineness Ratio 
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Height of Landing Gear 

ft 

g 

Gravitational Acceleration 

ft/sec 2 

K’ 

Inviscid Drag due to Lift 


K" 

Viscous Drag due to Lift 


K 

Drag due to Lift 


1 

Length of Airplane 

ft 


lh 2 


Liquid Hydrogen 
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Maximum Landing Weight 

lbs 

A N 

Subsonic Suction Parameter 
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Mean Aerodynamic Center 

ft 
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Location Of Aerodynamic Center 

ft 

r le 

Leading Edge Radious 
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Wetted Area of Wing 
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Exposed Planform Area 

ft 2 

S B 

Maximum Cross SEctional Area 
of Body 
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Thickness to Chord Ratio 


T.F.R.J 

Turbo Fan Ram Jet 


UK 

Kinetic Friction Coefficient 


V 

Landing Velocity 

ft/sec 

W 

Landing weight of Airplane 

lbs 

x cg 

Location of Center of Gravity 

ft 


Angle of Attack 

degrees 

& 

Glige Angle 

degrees 



Leading Edge Wing Sweep 
Trailing Edge Wing Sweep 


degrees 

degrees 


Taper Ratio 



Introduction 


The design group project chosen this year is to design a hypersonic 
aircraft which uses scramjets to accelerate from Mach 6 to Mach 10 and 
sustain that speed for two minutes. Initially the main difficulty was 
deciding on a propulsion system to propel the aircraft from the launched 
speed of Mach .8 to Mach 6 at 100,000 feet. The different propulsion 
systems being considered were solid rockets, liquid rockets, ramjets, and 


turoofan-rarrnets. The final decision was that the aircraft would use one 
full scale turoofan-ramjet. Later it was decided to add two solid rocket 
boosters to save fuel and help the aircraft pass through the transonic 


region. 

After the propulsion system was decided each member of the group was 
assigned a task such as: aerodynamics, aircraft design, stability and 
control, cooling systems, mission profile, and landing systems. The 
members researched their specific area of assignment and tried to get a 
physical understanding of how their system interacted with complete 
hypersonic aircraft. 


After two weeks of research the group was ready to begin to set the 
aircraft configuration. The two possible choices available for the 
configuration were a wavender or conventional design. The conventional 
design was chosen due to its landing characteristics and the relative 
expense compared to a waverider. It was apparent from the start that the 


amount of fuel required to complete the mission greatly effected the 
relative size of tr.e aircraft. Also the integration of the engine systems 
and their inlets affected the aircraft configuration. Each of the following 
sections of the report contribute to the final design of the hypersonic 
aircraft and is written by an individual member. 



GENERAL A iRCR A FT DES.IQM 


An unmanned hypersonic researched vehicle was designed to test a 
SCRAMJET. This vehicle had to meet the following performance 
specifications: 

1 ) The aircraft had to be carried to 40,000 ft. and launched from a 
carrier aircraft at M=.8. 

2) The aircraft will be capable of accelerating to M=6 and climbing to 
100,000 ft using its Internal engines. 

3) The aircraft will carry a prototype SCRAMJET engine that will be 
employed to accelerate the aircraft from M=6 to M=10 and is to 
maintain M= 10 for 2 minutes of study level flight. 

4) The aircraft will then decelerate and land back at its base, NASA 
Dryden Flight Test Center, CA. 

The aircraft design process started by choosing the propulsion system 
to meet performance specification 2. Solid rockets were considered 
because of their simplicity, high thrust, and the fact that they can be 
dropped to reduced drag, but they were found to weigh too much of the total 
vehicle weight. Liquid rockets were considered because of their high thrust 
but the fuel and oxidizer would require too much volume. Ramjets were 
considered because of their simplicity and light weight but were found to be 
lacking performance in the subsonic and high drag transonic region. Finally, 
a General Electric hydrogen fueled augmented turbofanramjet was chosen. 



The turbofanramjet was selected over the others for several reasons. It 
produces good thrust and efficiency within the Mach range of .8 to 6. It has 
an acceptable weight. It burns hydrogen, the same fuel as the scram jet. 
And In addition, the aircraft could have a powered landing. 

The turbofanramjet was selected full scale; 22.1 ft. long, 6.92 ft. in 
dia., and 6100 lbs. in weight. The reason for selecting it full scale was to 
assure enough thrust around M-1 since this region would have the highest 
drag and the least thrust available. 

Since the aircraft would burn liquid hydrogen, which has a very low 
density of 443 lb/ft' 5 , most of the aircraft volume would be taken up by 
liquid hydrogen. Thus, in order to get a preliminary sizing of the aircraft, a 
rough estimate of total fuel consumption was needed. This was done by 
approximating the average vehicle weight at 45,000 lbs. Next an average 
lift to drag ratio of 2.5 was assumed. From this an average drag of 
18,0001bs was calculated, Next a preliminary altitude versus Mach number 
flight path was selected for the turbofanramjet using the General Electric 
data. From this an average thrust of 60,000 lbs., an average fuel flow rate 
of 60,000 lbs, /hr., and a SFC of 1.0 was approximated. This resulted in an 
average accelerating force of 42,000 lbs. This then allowed the calculation 
of a flight time of 3 minutes from M~ 8 at 40,000 ft. to M=6 at 100,000 ft., 
which meant 3000 lbs. of liquid hydrogen would be burned by the 
turbofanramjet. From M=6 to M= 10 the SCRAMJETS would be turned on. An 
SFC of 2.0 and a thrust of 35,000 lbs. was assumed for the SCRAMJETS. The 
5CRAMJETS would then operate at M=10 for 2 minutes. From this the 
estimated SCRAMJET on time was calculated to be 7.48 minutes in which 
time they burned 7593 lbs. of liquid hydrogen. Thus the estimated total 
amount of hydrogen burned was 10593 lbs. This meant a fuel tank with 



2391 ft 3 was needed. 

A conventional configuration was next selected because the aircraft 
was hypersonic research vehicle for a SCRAMJET which would only fly for 
short amounts of time covering a Mach range from 0 to 10 and would not be 
at any one Mach number for a long period of time. A conventional design 
would also keep the plane simple and within current technology. 

A simple conventional design and the fact of wanting a slender low 
drag aircraft led to using one fuel tank placed in front of the 
turbofanramjet. Since some lift was wanted out of the fuel tank it was 
made slightly wider than it is high and is curved longer on the bottom. As a 
result the fuel tank was made 8 ft. wide , 6,5 ft, high , and 60 ft. long. This 
resulted in the fuel tank having 2450 ft. 3 which met the fuel requirement. 
The plane, however, resulted in being 90 ft. long and resembled a missile 
with a SCRAMJET strap to the bottom (figure D1.). 

After intial exact calculations were done on the plane in figure D1. it 
was found necessary to decrease the wing loading by one-half and add 
lifting surfaces to the front part of the aircraft in order to keep the nose up. 
This was done by increasing the wing span 10 ft. and strecthing the wing 40 
ft. up the side of the fuselage (figure D2.). The wings extending up the side 
of the fuselage were blended into the body like that of an SR-71. This also 
allowed for more fuel storage. The exact SCRAMJET data from General 
Electric also caused a reconfiguration of the SCRAMJETS' location. Since 
the SCRAMJETS produced small net thrust of on the average of 6500 lbs at 
100,000 ft. from M=6 to M= 1 0. It was decided to use two of them. And 
since they required large exit areas of on the average of 35 ft 2 they were 
placed on either side of the fuselage and under the wing so as to get some 
inlet compression from under the wing. The exit nozzles expanded under and 



over the wing and then expanded into the turbofanramjet nozzle to satisfy 
the necessary exit area requirements. The turbofanramjet nozzle was now 
made two dimensional and resembles a 6x7x10 rectangle. This was done 
in order to expand the SCRAMJET nozzles into the turbofanramjet nozzle. 
The full expanded SCRAMJET nozzle areas are shown in blue on figure D2. 
This expansion into the turbofanramjet nozzle also elimanted the base drag 
from the turbofanramjet nozzle when It was off. The base drag from the 
SCRAMJETS was eliminated by folding the nozzles into the side of the 
fuselage when not in use and expanding them outward when in use (figure 
D2.). The inlet to the turbofanramjet was placed on the bottom of the 
fuselage In order to assure a good compression surface for the inlet. The 
inlet to the turbofanramjet as well as the inlet to the SCRAMJETS will be 
closed off when not in use in order to reduce drag. 

After more exact calculations on the design in figure D2, the design in 
figure D3 was arrived at. It was found necessary to add two more 
SCRAMJETS in order to achieve an accelartion quick enough to keep the 
flight time and distance reasonable, ie., under 30 minutes and under 2000 
miles. The two SCRAMJETS were added by just making the existing 
SCRAMJETS twice as wide (figure 03.), The SCRAMJETS expansion nozzle 
area is shown in blue on figure 03. A detailed view of a SCRAMJET itself in 
both the on and off stages Is given In figure D4, The aircraft wings were 
also smoothed out to reduce drag (figure D3). The landing gear was also 
reposition in order to reduce flexure stress on the ground. The front nose 
wheel is keeped cool by the fuel tank. The avoinics is positioned in the 
front of the plane, 

This design (figure 03) resulted in aircraft 90 ft. long, 15 ft. -high, 
fuselage width of 8 ft, and height of 6.5 ft. It has a wing that has a span of 



40 ft., thickness to chord ratio of ,05, an area of 850 ft. 2 , aspect ratio of 
1.88, and a wingloading of 42.68 lb/ft 2 empty and 56.20 lb/ft 2 full. The 
aircraft has a gross weight of 47,774 lbs. and a empty weight of 36,274 lbs. 

The aircraft in figure D3, however, was larger than the French's carrier 
plane's limit by 5 feet, that is, their maximum plane length limit was 85 
feet. This meant that this design had to be reduced by 5 feet. This could 
have been brought about in 2 ways. 

The first way was to scale down the entire aircraft as it is now. At 
first a 50% scale down was studied. This evalution discovered some 
important characteristics about drag, lift, weight, and volume in scaling 
down the aircraft. Since when the plane is scaled down the surface area and 
cross sectional area go down by the square. This reduces the drag by the 
square. The inlet and nozzle areas are reduced by the square. This reduces 
the mass flow rate by the square and thus the thrust is reduced by the 
square, Since the drag and thrust are both going done by the square it seems 
relatively simply to scale down the aircraft. The catch is that the volume 
is going down by the cube. This means the fuel volume is being reduced to a 
greater degree than the drag and thrust. This means that aircraft will run 
out of fuel before completing its mission. The good part about the volume 
going done by the cube, however, is that the structural weight also goes 
down by the cube. 

The second way was to simply remove 5 feet of length from the 
aircraft's fuselage and readjust the plane accordingly. 



AIRCRAFT DESIGN II 


Gerry Fuerst 

The final aircraft configuration was the result of a long evolution of 
aircraft designs. Originally, the plane resembled a missile that was 
powered by a turbofan-ramjet (TFRJ) and a scramjet. However, the lenght 
had to be shortened and the wing area increased. The nozzles of the 
scramjets had to be reworked inorder to make them workable. The final 
design turned out to be that of a delta winged, tailess airplain (see figure 
Dll 1). 

The overall volume of the fuselage was dictated by the number of 
engines and the enormous fuel requirements. The plane uses one full-sized 
TFRJ to take it from Mach .8 to Mach 6. After Mach 6, the plane will use 
four full-sized scramjets to take it to Mach 10. Because of the large 
amount of fuel that is required to power these engines, the length of the 
plane was set at 85 feet (the maximum length that the French will allow). 

The scramjets are located at the rear of the aircraft- two on each side 
of the TFRJ. When the 2-D nozzle of the TFRJ closes, it provides two 
verticle expansion surfaces for the scramjets (see figure Dll 2). These 
expansion surfaces, along with pressure boundaries, will be used as the 
nozzles for the scramjets. 

The aircraft has three seperate inlets. A variable geometry inlet for 
the TFRJ is located on the bottom of the fuselage, and a mixed 
compression inlet is located on each side of the plane. These two mixed 
compression inlets supply air to the four scramjets. The scramjet inlets 
actually start about 18.5 feet upstream of the actual scramjet units. This 
is where the sides of the fuselage turn outward at a 6.5 degree angle. This 
turn in the fuselage produces an oblique shock that rests on the outer lip 
of the scramjets at Mach 10. This shock creates the initial compression 
of the air entering the scramjets. 



The planform of the aircraft is basically a delta shape. Because of the 
fact that the aircraft will fly at Mach 10, the wing was designed with a 
high sweep angle and a low aspect ratio. The surface area of the wing is 
relatively small. However, this area could be even smaller for a plane 
flying at hypersonic speeds, but due to the fact that the plane must be 
landed, it was not reduced any further. The following parameters apply to 
the planform: 

Wing Span, b = 40 ft 
Root Chord, Cr = 40 ft 
Tip Chord, Ct = 5 ft 

Mean Aerodynamic Chord, MAC = 27.037 ft 
Taper Ratio, A = 0.125 
Sweep Angle of Leading Edge,A l£ = 69 deg 
Sweep Angle of Trailing Edge,A re = 23 deg 
Exposed Surface Area, Se = 722.5 ft2 
Aspect Ratio, A = 2.215 
Maximum Thickness Ratio, t/c = 0.0370 
Zero Lift Angle of Attack, = o deg 
Aerodynamic Center, Xac = ,5(MAC) 

The following parameters apply to the vertical stabilizer: 

Surfac Area, Sv = 120.75 ft2 
Root Chord, Crv = 27.5 ft 
Tip Chord, Ctv = 7 ft 

Sweep Angle of Leading Edge,A t£v = 73 deg 
Sweep Angle of Trailing Edge,A re = 20 deg 




Xvs = 17.69 ft 
Zvs = 6.435 ft 



Weight Analysis 
Thomas Greiner 

To better estimate the hypersonic aircraft weight for all of its flight 
systems an analysis was performed using a computer program developed 
by both the United States Air Force and the National Aeronautics and Space 
Administration . The initial program, P.D.W.A.P (Preliminary Design and 
Weights Analysis Program) requested basic preliminary design data such 
as fuel weight, number of crew men, type and number of engines to act as 
an initial data file for a larger more comprehensive weight program. The 
larger program called W.A.A.T.S, (Weight Analysis of Advanced 
Transportation Systems) bases its findings on a data base containing 
information on aircraft already in existence. Not only does the program 
provide weight estimates for the total aircraft and its subsystems but it 
also sizes the aircraft for fuel and flight regime. 

When the program was run for the hypersonic aircraft, one problem that 
arose was that it did not allow for scramjet engines. Since there are four 
scramjets on the current design this could introduce substantial error, 
however given the engines provided by the program, the most reasonable 
choice was to describe the scramjets as ramjets. When the 
turbofan-ramjet onboard was entered into the program a highly 
exaggerated weight estimate for the inlet (over 200,000 pounds) was 
returned. This error was compensated for by considering the engine to be 
a ramjet and then adding to the weight of its inlet to adjust it to the 
actual weight of the turbofan-ramjet. 

The data returned by the program indicated that the initial weight 
estimates were very reasonable. The take off weight computed by 
W.A.A.T.S was 44,133 pounds, see figure weight 1 while the group used 
43,000 pounds , see figure weight2. That eight percent difference is 
acceptable, however the computer program underestimated for the thermal 
protection sytems so the higher number will still be used. The computed 



height of the aircraft matched the 8 feet that is used in the current design 
while the calculated span is 32 feet while the aircraft's current span is 
40 feet. 

Overall the W.A.A.T.5. program computed a very similar weight 
estimate for the aircraft and its subsystems reaffirming the intial 
estimates used by the design team. 



Rocket Boosters 
Thomas Greiner 

Since the main driving force in sizing the hypersonic aircraft was the 
need to hold large volumes of liquid hydrogen, solid rocket boosters were 
considered for the main propulsion system that would be used to 
accelerate the aircraft from the launch speed of Mach .8 to the Scramjet 
operating speed of Mach 10. The solid rocket was considered because of 
the high thrust and accelerations which would not be a concern to the 
unmanned aircraft, also the booster could be jettisoned after its use. 
Preliminary calculations indicated that the rocket weight would be 
75,000 pounds which was almost twice the aircraft weight. The idea of 
using rockets as the main propulsion system proved to heavy. 

A full scale turbofan-ramjet, (T.F.R.J.) was chosen to accelerate the 
aircraft fom Mach .8 to Mach 6 but a problem arose because the engine's 
burning of liquid hydrogen increased the aircraft's size. 

Another problem that occured was the efficency of the T.F.R.J. inlet to 
produce a reasonable pressure recovery at the low Mach numbers. 
Estimates made by the U.S.A.F. showed that a loss in thrust of 30% would 
occurr through the transonic region, where the aircraft drag increases 
sharply. There is some concern whether the engines could produce enough 
thrust to push the aircraft through the sound barrier and that it would take 
more fuel to do so because of the lower accelerations. 


These problems were great concerns that needed to be corrected 
without penalizing the design by increasing its size to add more liquid 


hydrogen. One solution was to put the aircraft into a dive and use gravity 


to gain speed after seperation from the carrier aircraft. This method 
could work but it 'would take more fuel to climb back to the original 
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Another solution was to use strap oq booster packs to accelerate 
through the transonic region. These boosters would be used from Mach .8 

to Mach 2.5 to accelerate the aircraft while the T.F.R.J. would produce 
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enough thrust to overcome drag. This T.F.R.J. was not used at full thrust to 
reduce the amount of fuel needed to span the transonic region but it was 
kept operating to allow for a smooth propulsion transition when the solid 
boosters were jettisoned. 

The two boosters are integrated into the aircraft body, one under each 
wing, see figure Booster 1 mounted next to the scramjets. The booster's 
nose is slanted to lower the drag and also the two-dimensional exit nozzle 
is located partially behind the aircraft body to reduce the drag. The length 
of each booster is 15 feet and their combined weight is 15,000 lb. This 
additional weight puts the aircraft's takeoff weight at approximately 
62,000 pounds. 

The size and weight of the booster rockets was calculated using a 
spreadsheet, see figure Booster2 which was based on the constant 
acceleration that was assigned. An acceleration of 1-6, (32.2 ft/s'2) was 
chosen to reduce the range traveled and the amount of fuel used. Using 
the given acceleration the amount of fuel required to maintain the chosen 
flight profile is calculated and then the required volume is determined. 
The solid propellant chosen was J.P.N. which has a specific impulse of 240 
seconds. This was the upper limit for a solid rocket fuel impulse. The 
initial -weight of the aircraft was 62,000 pounds and the weight of the 
fuel burned by both the solid rocket boosters and the T.F.R.J. were 
subtracted from the aircraft weight as it accelerated. The thrust 
required to accelerate the adjusted mass was calculated and then 
converted into pounds of propellant. Knowing the density of the 
propellant the volume required to contain the propellant was determined. 
With this information a rocket diameter of 2.5 feet was determined to be 
exceptable in terms of drag and volume for the rocket's tube. 

The overall benefits of using the booster pack is a reduction in range of 
: miles due to the large acceleration, a reduction of 500 pounds of fuel 
needed for the T.F.R.J. in the transonic region, and the additional thrust to 
offset the thrust lost due to an inefficent T.F.R.J. subsonic inlet. 



CENTER OF GRAVITY LOCATION 


Gerry Fuerst 

The center of gravity for the entire aircraft was found by taking the 
weighted average of the center of gravity for individual aircraft 
components. The Analysis was as follows: 


COMPONENT 

WEIGHT (LBS) 

LOCAL C.G. (FT) 

Nose * 

3365 

77.5 

Body 

5342 

40 

Fuel Tank 

13913 

49 

TFRJ 

6500 

20 

Scramjets 

12000 

20 

Inlet 

400 

40 

Front Landing Gear 

50 

50 

Rear Landing Gear 

150 

1 1 

Wing 

4675 

23.125 

Vertical Tail 

664.125 

1 1.14 


v Mote: the nose section contains a portion of the main fuel tank. 


V/hen the fuel tank is empty, it not only affects the weight of the fuel 
tank but also the weight of the nose. Therefore, the empty weight of the 
fuel tank and nose is 4413 lbs and 1365 lbs respectively. With the above 
information, the center of gravity was calculated as follows: 



Sum (Weight*Local C.G.) 


Xcg = 

Sum of the Weight 

Xcg (full tank) = 35. 102 ft 
Xcg (empty tank) = 29.053 ft 


The C.G. is measured with respect to the rear of the aircraft. 



nignc nrome i- 
Thomas Greiner 

A necessary requiremet for any hypersonic vehicle is to minimize the 
size and weight of the aircraft. The size and weight requirements for this 
mission were specified by the French based on their launching aircraft 
requirements. Since the hydrogen fuel used has a large storage volume it 
is necessary to minimize the amount of fuel required for the flight profile. 
This mission profile consisting of a five phase flight was initially used to 
minimize fuel. 

Phase One) Acceleration from Mach .8 to Mach 2.5 at 40,000 feet. 

Phase Two) Following the Q-Curve of 1850 at Mach 2.5 at 40,000 feet 
and accelerating to Mach 6 at 75,000 feet while 
maintaining a climb angle of four degrees. 

Phase Three) Leaving the Q-Curve at 75,000 feet and climbing to 
1 00,000 feet whi le maintaining Mach 6. 

Phase Four) Leveling off at 100,000 feet turning on the scramjets, 
turning off the turbofan-ramjet and, accelerating to 
Mach 10 and sustaining that speed for two minutes. 

Prase Five) Turning off the scramjets and slowing down below Mach! 

and turning the turbofan -ramjet on and heading to Dryden 
Air Force Base m California under power. 

The preceding phases were analyzed using a spreaaheet, see figure 
Profile! which uses conditions such as acceleration, time-of -flight, and 
the climb angle to minimize fuel consumption. Calculations were 
conducted to follow a constant dynamic pressure line, a O-curve, which 
will allow the aircraft to maintain a constant angle of attack, from 40,000 
to 100,000 feet. Originally, a Q based on Mach 6 at 100,000 feet was 
attempted Put could not be used due to two factors, the magnitude, of the 
aerodynamic fo r ces at the lower altitudes and the lack of thrust of the 
turocfar.-ramjet at higher altitudes to meet the required accelerations 
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altituded required to stay on the Q-curve for the entire flight. The Q value 
is above 2000 which implies there would have to be a substantial increase 
in structural weight to accept the loads. For a hypersonic aircraft a Q 
should be about 1 800 therefore a new value was chosen based on Mach 6 at 
75,000 feet with a Q of 1850. This an acceptable value that will not 
greatly effect the intial estimate of the structural weight. 

The aircraft follows this Q-curve in phase two until it reaches Mach 6 
at 75,000 feet. At this time the aircraft enters phase three were it 
leaves the G-curve while maintaining a climb angle of four degrees and 
climbs 100,000 feet while sustaining Mach 6. In this high altitude period 
of the flight profile the thrust of the turbofan-ramjet falls off, however 
the thrust required also drops. Since there are no other acceleration 
requirements to meet other than those caused by the change in the speed 
o"' sound due to the atmospheric temperature, at this point the thrust has 
only to overcome the drag and the force from the weight component, 
therefore the thrust available is still large enough to meet the thrust 
required. 

in phase four of the flight the aircraft levels off at 100,000 feet, the 
turbofan-ramjet is turned off and the scramjets are turned on. Based on 
the data supplied by the spreadsheet it became obvious that the aircraft 
needed more thrust to reduce time of burn and, in turn, the fuel required to 
accelerate from Mach 6 to Mach 10. Two additional scrarnjet modules 
were added to the aircraft. Adding the engines increased the thrust and 
allowed for considerable increases in acceleration and reduced the amount 
of fuel burned during this phase of the flight. 

Finally in phase five after the two minute engine test the scramjets 
a-'e so ut c'f ana the plane is decelerated by the force of drag. The aircraft 
slows down to Mach 1 at 50,000 feet and the turbofan-ramjet is turned on 
while the aircraft continues to slow down. The hypersonic aircraft 
maintains Mach .8 and flies to Dryden Air Force Base in California and 


lands. The actual flight path along the west coast of the United States is 
illustrated in figure Profile2. 

To further reduce the liquid hydrogen necessary to accomplish the 
mission an energy state program, was written. 



MINIMUM FUEL FLIGHT PATH 
(ENERGY STATE METHOD) 
by 

Bob Stonebraker 

The accelerated accent from release to Mach 10 at 100,000 feet 
was broken into two phases based on the engines used. Phase 1 
employs the turbofan-ramjet and two booster rockets to reach Mach 
6 at approximately 76,000 feet (q = 1800 psi.). Phase 2 then 
continues the accent under power of the scramjets to Mach 10 at 
100,000 feet to begin the two minute test. The energy state method 
showed that the altitude vs Mach number flight path which consumed 
the least fuel was simply one of a constant high dynamic pressure. 
Temperature limitations determined the highest dynamic pressure to 
be 1800 psi. 

The data provided for the General Electric augmented turbofan- 
ramjet engine consisted of net thrust (Fn) and specific fuel 
consumption (sfc) for various Mach numbers and altitudes. To make 
this data easier to incorporate into a computer program, it was 
used in a subroutine (SUBROUTINE TFRMJT) to produce values for any 
Mach number at any altitude. To do this, Fn and sfc were 

separately plotted vs Mach number for a constant altitude. This 
was repeated for various altitudes yielding several curves, one for 
every 10,000 feet from 40,000 to 100,000 (Figure ES-1) . A plotting 
package was then used to obtain best-fit equations of these curves 
which modeled the given data quite accurately. These polynomials 
were then used in the subroutine in conjunction with the cubic 
spline interpolation method for values at altitudes between those 
of the equations. In this way, the subroutine yields a 
corresponding net thrust and sfc for any Mach number (.8-6) at any 


\i 



altitude (40,000* - 100,000'). Where several values of Fn and sfc 
were given for the same Mach number and altitude in the original 
data, that for the best fuel economy was chosen for the plots. 

The given data for the experimental SCRAMJET was also reduced 
in a similar manner through a subroutine (SUBROUTINE SCRAMJET) 
which yields values for a range of Mach numbers (6-10) at any 
altitude (80,000' - 100,000'), (See Figures ES-2 & ES-3) . The data 
provided, gave thrust (Fn) and specific impulse (Isp) values as a 
function of dynamic pressure and fuel/air equivalency ratio for 
various values of altitude and Mach number. Because the data for 
this engine was given in terms of various dynamic pressures (q) 
instead of altitude, the best-fit equations in this subroutine have 
altitude as the independent variable with Mach number held 
constant. Then the cubic spline interpolation subroutine 
(SUBROUTINE SPLINE) finds the desired net thrust and Isp for non- 
integer Mach numbers from 6 to 10. Once the Isp (sec.) is found, 
the sfc is computed by the following; 


SFC - ^°° 


Isp 

Two additional subroutines had to be written for the Energy 
State program; one for standard atmospheric data as a function of 
geometric altitude (SUBROUTINE ATMOS FR) and another for the total 
drag of the aircraft as a function of Mach number and altitude 
(SUBROUTINE DRAG) . The total drag data, computed by another team 
member, was received in tabulated form and was therefore also 
reduced through curve-fit polynomials and cubic spline 
interpolation. 



The energy state method can be used to determine the minimum 
time or minimum fuel required to reach a Mach number and altitude. 
To compute the minimum fuel accent required five separate programs. 
The first code (PROGRAM HECONST) determined lines of constant He 
(ft) , which defines the amount of potential and kinetic energy that 
an aircraft possesses at a certain Mach number and altitude. He 
is given by; 






Two additional codes determined contours of constant Fs for phase 
1 and phase 2 of the mission (PROGRAM FSCONST & PROGRAM FSCNST2) . 
Fs is the vertical distance traveled per pound of fuel burned and 
is given by; 

p 5 pA ol(Tn~ b) 

Fs ' xTsfc' TVS^c- w" 

The minimum fuel trajectory is one in which the lines of constant 
Fs are tangent to lines of constant He. These plots showed that 
a greater thrust per pound of fuel is obtained at high Mach numbers 
and low altitudes. However, aerodynamic heating obviously limits 
flight in this realm. Therefore, the minimum fuel accent was 

chosen to be one of a constant q = 1800 (Figure ES-4) . For 
comparison, the quantity of fuel required to accelerate from Mach=6 
to Mach=10 all at 100,000 feet was determined to be approximately 
25,000 lbs. This is obviously unacceptable. 

Once the flight path was determined, several values could be 



computed along the trajectory. Figure ES-5 shows the thrust 
available and required curves. The last two codes (PROGRAM PHASE1 
and PROGRAM PHASE2 ) were used for this and to compute the fuel 
consumed and elapsed time for phase 1 and phase 2 of the mission. 
The fuel burned is given by; 

AkJy - j — J. 

The elapsed time is given by? 



The range required for these accents was also computed, as shown 
below, and found to be approximately 130 miles for phase 1 and 970 
miles for phase 2. 

4 - J(tf a. iff - (a 

Originally, the plan was for Edwards A.F.B. to be the base of 
operation over all phases of the mission, from carrier aircraft 
takeoff to test vehicle landing. The plan incorporated a 
triangular pattern with one leg over the ocean for the accent, 
acceleration and scram jet test phases. The aircraft was to be 
carried to the release point approximately 300 miles southwest of 
base. This would allow a return, either to Edwards or Vandenburg 
A.F.B. under power of the turbofan-ramjet in the event of scramjet 



misfire. However, after a better estimate of required acceleration 
distances, it was necessary to alter the plan. It was decided to 
originate the mission from Seattle and release the aircraft over 
the ocean, then run approximately 100 miles off the coast to land 
at Edwards. This plan still fulfills the requirement of non- 
supersonic test vehicle flight over populated areas. 

Control of the mission from release to landing is to be fully 
automated. Based on this, it was determined that the avionics 
should include a receiver, transmitter, inertial navigation 
equipment, stability augmentation systems, and a computer. This 
equipment was estimated to weigh approximately 1000 lbs. The 
aircraft cooling system was not included in the avionics. The 
control program will originate from systems based at Edwards and 
be communicated to the aircraft via satellite. Real time 
communication between the aircraft and Edwards A.F.B. is also 
necessary to monitor data from the test scramjets. The on-board 
computer will maintain a backup control program in addition to its' 
stability augmentation and engine monitoring functions. 
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AERODYNA MICS 


The lift and drag of an airplane are the key factors In determining the 
type of engines and the performance of the airplane. For the first phase of 
the design process , an initial configuration was approved, an estimate of 
the CD 0 of our airplane was calculated using the build up method. The 


initial calculations indicated that using one full scale GE turbofan-ramjet 
would be sufficient to overcome the drag of the airplane at the various 
phases of flight where GE's turbofan-ramjet would be used. - After ther first 
phase of the design process was completed, the drag of the airplane was 
more accurately calculated using Nicolai’s (chapters 2 & 1 1 ), and the 
Datcom reference. After several iterations on the drag numbers because of 
slight modifications to the airplane design , final drag numbers were 
calculated to the final airplane configuration. This data was then given to 
the flight path group which was responsible for determining the flight path 
and the fuel requirements of the mission. After the fuel requirement and 


time of flight for the mission profile were calculated using the thrust 
provided by the two scram-jet engines it was obvious that an additional 
two scrarn-jets, one on each side of the airplane, would be required to cut 


the fuel requirement and time of flight 


These additions effects on the draq 


were included as a \ 0 % increase to the CD n of the airplane because these 


modifications were done the day before the final presentation for Winter- 
quarter. Therefore, the effects of these additions had to be calculated more 
a c c u r a t e i y i n 5 p r 1 n q q u arte r . 



The last minute modifications, that were introduced at the end of 
Winter quarter, were studied more closely at the beginig of Spring quarter 
and it was obvious that some modifications to the shape of the airplane had 
to be done. After the final shape of the airplane was modified and approved, 
the final drag numbers were calculated . The drag numbers for the old 
configuration could not be modified to account for the changes in the 
airplane configuration mainly because the fundemental parameters used in 
calculating the old drag numbers were drastically changed by the new 
design; these parameters include the airplane fineness ratio (lg/d), the 

nose fineness ratio (l N /d), the after-body fineness ratio (l A /d), the aspect 


ratio of the wing, the surface area of the wing, the taper ratio, and finally 
the length of the airplane, Therefore, the drag numbers for the new 


configuration had to be calculated from scratch. 


'he three coefficients 


f K -n #■ 
U kit 


had to be calculated to provide us with sufficiently accurate drag numbers 


were CD 0 , K (factor of drag due to lift), and CL . These three coefficients 


need to be calculated seperately for the four flight regimes in our flight 
path. 


SUBSONIC 


For the wing, the zero lift drag is composed of two parts, skin friction 
drag and pressure drag. As for skin friction drag, it is caused by shearing 


stresses within a thin layer of retarted air on the surface of the 
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the boundary layer. The amount of viscous resistance depends on whether 
the flow in laminar or turbulent; for our airplane turbulent flow was 
assumed for all flight regimes. As for the pressure drag, it is usually small 
compared to skin friction and it is primarily caused by the displacement 
thickness of of the boundary layer. Methods for predicting subsonic CD 0 for 

wings are essentially empirical and are based on streamwise airfoil 
thickness ratio (t/c). 

For the body, at subsonic speeds the drag of smooth slender bodies is 
primarily skin friction. The Reynold's number is based on body length, 
boundary layer condition, and surface roughness. The pressure drag is also 
generally small for fineness ratios above 4( The airplane's finess ratio is 
7.083) but becomes significant for blunt bodies. At the subsonic reigon, A 
CD 0 of 0.0129885 was calculated at M=0.8. Also at M=0.8 and at 40,000 ft, 

the drag is 7508.25 lbs with a L/D of 6.433. 

TRANSONIC 


For the wing, the transonic range varies greatly with airfoil shape and 
thickness, but for simplicity it can be considered to begin at approximately 
M=0.9 and end at M= 1 .2. Because of the mixed flows, drag in the transonic 


reigon does not lend itself to theoratical or experimental 


analysis. 


! 


wing drag at the transonic reigon is mainly composed of skin friction and 
wave drag. As for skin friction, a little increase in drag is experienced due 
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to viscocity . Therefore, the skin friction drag wil! be assumed constant 
and equal to the subsonic skin drag throughout the transonic range, As for 
the wave drag, the variables involved In a wing design that effect the 
manner in which shock waves develop on the surface are many; they include 
sweep, aspect ratio, taper ratio, thickness ratio variations between root 
and tip thickness, position of the maximum thickness, incidence and leading 
edge thickness. As a result its very hard to predict the wave drag in the 
transonic reigon. Estimates for the wave drag can be obtained by finding the 
drag divergence Mach number and then using figure 1 1 . 1 0 in Nicolai's to find 
CD 0 . 

For the body, the genera! approach consists of predicting the skin 
friction, the drag divergence Mach number, the variation of base drag with 
Mach number, and the variation of pressure drag for Mach numbers above l. 
Therefore, the drag of the body at transonic speeds consists of skin friction, 
base pressure, subsonic pressure drag, and supersonic wave drag, At the 
transonic reigon, as predicted, a CD 0 of 0.039493 at M=1,2 was the largest. 
Also at M= 1 .2 at 40,000 ft, the drag was 146983.92 lbs with a L/D of 3.273. 


SUPERSONIC 


For the wing, at supersonic speeds an increase in the Mach number 
results in a decrease in the skin friction coefficient at constant Reynolds 
numbers. This variation is primarily due to the variation in the temperature 



and density at the surface. The full reduction in skin friction at supersonic 
Mach numbers is justified only when stabilized conditions and zero heat 


transfer are attained: For transient flight the skin friction will be assumed 


equal to the incompressible value although in reality it varies between this 
value and the zero heat transfer. Another important factor for estimating 
drag at the supersonic speeds is the Reynolds number, this can be accounted 
for by taking into account the ratio of compressible to incompressible skin 
friction coefficients. As for the wave drag, the well known linearized 
supersonic theory is used in predicting wing wave drag. For the airplane the 
wave coefficients were obtained using equation in Nicolai's which are based 
on the supersonic linear theory. 

For the body, the characteristics of compressible skin friction drag for 
bodies are similar to those of wings; the skin friction coefficients decrease 
as the Mach number increases. As for the wave drag, two methods for 
estimating the fore-body and after-body wave drag are presented in the 
Datcom reference. The second method, which was used, is based on 
similarity parameters. The wave wave drag is seperated into the fore-body 
drag, the isolated after body drag, and the interference drag of the 
fore-body and center section of the after-body. In the supersonic reigon, a 
CD 0 of 0.0242 was calculated at M=3.0. Also at M=3.0 and at 65,5000 ft, the 
drag was 40666.36 lbs with a L/D of 1 . 1 592. 



HYPERSONIC 


For the wing, due to the non-linearity of hypersonic flow, approximate 
methods for estimating force characteristics are very desirable. Among 
the methods used, Newtonian and modified non-Newtonian flow theory have 
proved very useful. Newtonian theory is based on the assumption that the 
shock concides with the wing surface and no friction exists between the 
wing and the boundary layer. The fluid particles ahead of the wing are not 
disturbed until they encounter the wing. 

For the body at hypersonic speeds the drag of the body is caused 
primarily by the pressure and skin friction drag of the nose. Both the 
after-body and the base drag become insignificant at higher Mach numbers. 
The drag coefficient was calculated at M=6.0 and at M=10.0 using Newtonian 
flow theory and it was also calculated using the supersonic linear theory; 
the discrepancy between the two methods was less than 10% (the Newtonian 
flow provided the lower drag values). It was deccided to stick with the 
CD 0 s found using the supersonic linear equations since the discrepancy 


between the two methods was relatively small. At M=6.0 a CD 0 
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'was calculated and at M=!0.0 a CD 0 of 0.01551 7 was calculated. At M=6.0 


and at 76,000 ft the drag was 32294 lbs with a L/D of 1.4002, and at M=I0.0 
at 100,000 ft the drag was 26036 i lbs with a L/D of 1.449. 



After the CD 0 was calculated for the various flight regimes a 10% 

increase was added to CD 0 to include the interfernce effects. In addition, 

the coefficient of drag due to lift (K) had to calculated in order to account 
for the effect of lift on the drag of the airplane. For the subsonic reigon the 
equations in Nicolai's (Chapter 1 1) were used. For the supersonic and 
hupersonic reigons the K values were calculated uning the supersonic linear 
theory. Finally the values of K for the transonic reigon had to be 
approximated. After obtaining the CD 0 's and the K values for the various 

Mach numbers, the drag at the various altitudes and Mach numbers was 
calculated. Three seperate programs were written to calculate the drag and 
angles of attack at the various Mach numbers and altitudes. A listing of the 
equations and procedures used in calculating the drag are listed in Appendix 
A and the listing of the programs used in calculating the drag can be found 
in Appendix B. The tables AD1-AD6 and the figures ADI -ADS include the 
results of this chapter. 



BASIC STABILITY ANALYSIS 


Gerry Fuerst 

Due to the large amount of fuel on board the test plane, there will be a 
considerable shift in the C.G. during flight. The C.6. will shift down and 
towards the rear of the aircraft as the fuel is burned. The distance 
towards the rear of the plane that the C.G. will travel was computed as 
follows: 


C.G. Travel = 35. 1 02 - 29.053 » 6.05 ft 

For stability, it is necessary to have a positive static margin. The SM 
is directly related to the distance between the center of gravity and the 
aerodynamic center. Although the C.G. travels back 6.05 feet, it always 
remains ahead of the aerodynamic center. This means that the static 
margin always remains negaive. The calculation of the static margin was 
as follows: 


Static Margin: SM - (Xac - Xcg)/MAC 

Full Tank: SM = 13/27.037 = 0.4808 
Empty Tank: SM = 6.875/. 037 = 0.2543 

For an aircraft to be statically stable, its value of Cm(alpha) must be 
negative. When this is the case, the aircraft is trimmed at a positive 
angle of attack. When the angle of attack is suddenly increased, the 
aircraft will generate a negative moment to push the nose back down to 
the origional trimmed angle of attack. The Cm(alpha) for the test plane 
was found to be negative for both the full and empty tank cases. It was 

r- r% r-r-i } » f r y r\ -or- f t\ 1 1 » i > r > , 
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Cm(alpha) = -CL(alpha)*(Xac - XcgVMAC = -CL(alpha)*SM 



Case 1 : Mach .8 


Full Tank: Cm(alpha) = -.7965/rad 
Empty Tank: Cm(alpha) = -.4213/rad 

Case 2: Mach 10 

Full Tank: Cm(alpha) = -.1957/rad 
Empty Tank: Cm(alpha) « -.1035/rad 

The Cm eg of the test plane was calculated for three different cases: 
subsonic (M=.8); supersonic (M=2); and hypersonic (M=10). Refer to tables 
SI thru S3 for Cm eg versus angle of attack, CL, and CD. Cm eg is plotted 
versus angle of attack in figures SI thru S3. 
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Gerry Fuerst 


The landing gear of the test plane will consist of two rear skids and 
one forward wheel. This configuration was choosen pr-imarilly for its 
reliability and simplicity. The rear skids are better suited in handling the 
high temperatures that will be experienced during flight. A forward wheel 
will be used because it will allow for the steering of the aircraft after it 
has landed. This wheel will be kept cool in flight by surrounding its 
compartment with the liquid hydrogen fuel tank. The rear skids are 
located 1 1 feet from the rear of the aircraft, while the front wheel is 50 
feet from the rear. 

In designing the landing gear, it is necessary to determine the 
maximum loads that will be experienced during landing. These maximum 
loads will occur as a short impulse when the landing gear touches down. It 
is assumed that the test plane will land at a velocity of 21 1 ft/s, and it 
will be descending with a glide angle of 5 degrees. The glide angle is the 
angle that the plane's descent makes with the ground. The 5 degree glide 
angle will be the maximum that the plane experiences and was chosen to 
determine the largest possible loads on the landing gear. Finally, the last 
assumption that was made is that the time of impulse will be 0.5 seconds. 

A simple free-body diagram of the landing gear is as follows: 



f 
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Constants: 


Landing Weight of Plane, W - 47000 lbs 
Height of the Landing Gear, h - 4.0 ft 
Kinetic Friction Coefficient, Uk =0.6 
Landing Velocity, V = 21 1.0 ft/s 
Glide Angle, Theta = 5 deg 
Time of Impulse, dT = 0.5 sec 
Gravitational Acceleration, g = 32.2 ft/s2 
The maximum force, on both skids, in the vertical direction (Fy) can be 
calculated from the following formula: 

(W/g)*Vy 


Py = where Vy = V sin (Theta) 

dT Vy- 211.0 sin (5) 

= 18.39 ft/s 


(47000/32.2)* 1 8.39 

py- 

0.5 


Fy = 53685.093 lbs 


This is the maximum vertical load that will be experienced by both 
skids together. Each skid will, therefore, experience half of Fy or 
26842.547 lbs. 

The maximum force experienced by both skids in the horizontal 
direction (Fx) can be calculated from the following formula: 
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Fx - Uk * N 
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where N is the maximum normal force 

The maximum normal force on both skids is equal to 

the maximum vertical force on the skids (53685.093 lbs). 


Therefore-, 

Fx = 0.6 * 179279.503 lbs 
= 3221 1.056 lbs 

The maximum horizontal force on each skid seperately will be half of 
Fx or 16105.528 lbs. . 

The maximum moment produced by each skid (Mz) can be computed from 
the following formula: 

Mz = (Fx/2) * h 

= (16105.528) * 4.0 ft 
»64422.112 ft lb 


The following is a free-body diagram of a skid with the maximum loads 
in place: 
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Materials and Cooling Systems 


A A 


Throughout nistory, tne advancement of Higher speed aircraft has 
been plagued with tne problems of what materials, size and sfrencitn, anc 
cooling systems to use, that are most beneficial to tne mission High 
thermal loadings ana body forces acting along tne entire aircraft require 
new technologies to be developed This section of the report deals with 
ideas anc concepts mat can oe usee m the development c? this aircraft. 

The first step m developing tne materials ano cooling systems is to 
calculate a temperature distribution aiong tne aircraft, i r»e temperature 
distribution is snown through color cooing in Figure MSP The program ano 
data tables are snown m appendix D. This program is cased on empirical 
data and the geometric structure of the aircraft, it should be noted that 
the program created is for the worst case scenario Tms program does not 
include soanwise nearing, this is tne heating aiong or parallel to tne 
leading edge. The values in the program were found in many Doors anc 
converted to nave tne same units. The emissivity of tne material will De 
considered to be .89, the reason will be discussed later Note that in the 
Figure MSI it can oe seen that tne boundary layer cnanoes from laminar to 
turbulent at approximate 5.5 feet (noted m tne mqner temperature) Tnis 


program does not however consider the ‘transition 
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turbulent. ; ne hottest temperatures can be seen to oe at me ieaainQ 
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edges, the nose, the inlets, and the vertical tail . 
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these temperatures can oe cocaines Noting me mgr temperatures, me 
first materia! to Iook at would oe a ceramic because of ns hi tin resistance 
to temperature and neat. The problem ;s that ceramics a r e not fully 
developed. A ceramic material can withstand tne nigh temperatures, out 
the dynamic forces on me posy as great as they a r e. win carnage me 


ceramic, maybe causing mi lure cue to cracMnc 
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aircraft nave to have certain, characteristics. Tne material for rms 
aircraft snouio nave nigh thermal resistivity, low creed, and low oxidation 
coefficients Tne material cnosen is an alloy of aluminum and titanium. 
This maternal is very strong, but light, and has coos mams on the afore 
mentioned character isticsCftef 2.. Tne emissivity of this marena 5 r- ~ 
out with a special black coating, like on me 5R-7 1 , tne emissivity rises to 


.89. 


The material on me skip will averace in thickness from . i in. to .2 m 


The titanium ailumimde can withstand temperatures of up to 1500 F. 
which is very beneficiaURef. 1 X because as Figure m $3 shows, this 
substantially reduces me amount of wine area that neeos to ce active v 
cooled. 


Knowmq wnat types of materials to use a no wnere to actively ron* 


ORIGINAL PAGE IS 
OF POOR QUALITY 


AO 



the next objective is to pick, a cooling system to use, Narrowing down me 
choices them seems to be three useful choices (Fioures MS3.fiS4.MS5). 
Fiaure MS6 snows a panel corn' i duration for the inlet This is tne oniv area 

V V N 

along with tne unaer oeiiy of tne plane needed to oe cooiea oy 2 cane- 
configuration. The idea is that liquid nyarogen is pushed tnrougn tne 
panel, the nyarogen then acts as a neat sink, extracting neat from me skip, 
tnen the hydrogen is run tnrougn a neat exchanger where it is men let out 
into the free stream o r out mm the engine as fuel. This win only n* iisso 
tor the inlets ana unaer oeiiy of tne plane, it is not efficient enough for 
the entire aircraft Tne major advantage of this system is tne low weight 
its major disadvantage oeing mat it requires a large fuel neat sink. 

The nose is very complex when it comes to heating. Tne diameter of 
the nose and the curvature of tne cone all add to the complexity c f me 
system Therefore, tne use o T ' 3 nose cap is the proper choice ■mis rose 
cap is composed of a JTA carbon composite This cap is verv strong arc 
durable, plus me reusability of tms cap is a major advantage, cost wise 
The leading eages require the most work, and are tne major problem. 
The wines ana vertical tail require neattna with 3 larae system Tne 




eading edges can oe cooked by two methods. Figure MS 4 snows a sora* 


cooling system, me 00 1 active of ims system is that a Maine is to oe 


t ^ r*. 
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exchanger, the system then starts over again in a continuous cycle, me 


major. advantage of this type of system is that the leading edge radii can 
be very small, out the major disadvantage is that it is unforgiving to local 


failure. 

The second type of leading edge system is a tube system, Figure MS3. 
The concept behind this system is that the leading eage would be composed 
of a carbon-carbon composite which is meshed or embedded with tubes 
that run throughout the leading edges. The tubes filled with a liquid run 
through the leading edge acting as heat sink, then me liquid is transported 
to a heat exchanger where the cycle starts all over again. The advantage 
to this system is that the cooling system can act oast the leading edge, 
even to six feet oast the leading edge, onto the wings where it. is needed. 

These two types of cooling have good points, but the deciding factor 
is that in the caroon-carbon composite tube cooling tne idea of cooling 
past the leading edge outways any aspect that the spray cooling has. Tne 
carbon-carbon tube cooling will act for the first fifteen feet of the wing, 
after that, areas that needed cooling can be cooled without carbon-caroon, 
just tube cooling will be used af'ter(Ref. 3). 

The vertical tail is quite similar to the leading edges of the wing, 


except that the maxirnuirn temperatures expected will be on the order o 


2000 F, The tail will be cooled with lust a simple tune coonna svstern. 
with carbon-caroon lust at the first foot of me vertical tail 
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The fuel tanks are very difficult to keep at the proper temperature. 
On one side of the fuel tank the temperature is approximately -419 F, 
while on the otherside the temperature is approximately 1800 F. There 
are two types of purge systems that can be used. 


The first purge system is a N 2 purge system (Figure MS6), The major 

advantage of this type of system is that you can use available insulations 
and it is reasonably inexpensive, Its major disadvantages is that it has 


temperature limitations and it is difficult to inspect. 

The second purge system is the C0 2 purge system (Figure MS7). The 

major advantage of this type of system is that it limits liquid hydrogen 
boil off and there is no liquid phase. Its major disadvantages is that there 
are complex ground handling requirements and it to is difficult to inspect. 


'ne decision is to use the COo system basically because the carbon 


dioxide is better suited for the heating and cooling, and also because it. 
limits liquid nydrogen boil off; Liquid Hydrogen boil off is of major- 
importance in deciding what to use, because you could ultimately lose 30 
percent of your fuel to hydrogen boilofl The addition of a vapor shield on 
the tank wall will further decrease the boiloff inside the tank(Ref. 3). 

The choice of which coolant to use for the cooling systems for the 
leading edges was simple to choose The only liquid with a high specific 
heat so as to absorb more energy is lithium. The other liquids in 



comparsion are not even close, it will be easy to use and is the most 
efficient. 
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En el Ta nk in sulatio n Thickne ss-Fjm a 

Mach 10 Hypersonic Aircraft 
By John Oon T ong & Gerry Fuerst 


In hypersonic flight, it is necessarily to hove to hove sufficient 
insulation around the tank to maintain the hydrogen in a liquid form. A 
slight increase in temperature of liquid-hydrogen fuel may results in 
tremendous increase in pressure. For example, hydrogen at its liquid 
state at -4232F has a density of 70.1 g/liter but at its vapor state, it 
has 1.3 g/liter. If insulation is not used adequately, the tank will 
explode. It is our goal to design an insulation system that will 
insulate the fuel used for a cruising time of approximately 20 
minutes. 

The derivation of the equation used to calculate the thickness 
of insulation was taken from NASA Technical Memorandum Faoer Th 
X-2025 by Mark D. Ardema. The author used standard analytical 
techniques to develop procedures for calculating the insulation's 
thickness. The assumptions that were used are : 

1. Only heat transfer by conduction is considered since radiation and 
convection effects are negligible compared with conduction. 

2. Thickness of insulation used is much smaller than the diameter of 
the fuselage. 

3. Conductivity of insulation is much smaller than the conductivity of 
all structural elements. 



4. Thermal constants are independent of time, temperature, and 
position. 

5. Insulation is continuous and homogeneous. 

6. The temperature at the wall of the inner fuselage is the same as 
the temperature of the liquid-hydrogen fuel. 

7. Exterior surface of fuselage is exposed to square temperature 
pulse. 

With the above assumptions, the following equation is obtained: 

hf^/fcH - Ktf(Tg-T||)/L 2 "- 2K(TS-T«)/k *(1/12 + 2 (-1)" (l/n 2 * 2 + 
2kt f /L 2 ) ) = o 

where the definition of the symbols used are as follows: 
h^ - hydrogen heat of transformation, (Btu/ft*3) 
k - diffusivity, (ft*2/hr) 

K - insulation conductivity, (Btu/hr-ft-SF) 

Lb - tank thickness,(ft) 

tf - cruise time,(hr) 

T^ ~ liquid hydrogen fuel temperature,(2F) 

T 0 - initial exterior surface temperature, (2F) 

T s - cruise exterior surface temperature,(2F) 

® - insulation density,(lb/ft*3) 

® B - tank density, (lb/ft‘3) 



- hydrogen fuel dens1ty,(lb/ft*3); 

It can be seen from the equation that for fixed materials and 
temperatures, the optimum insulation thickness is proportional to the 
square root of the cruise time. The above equation can be simplified 
further if we consider the steady-state condition. The steady-state 
equation can be obtained by setting the initial exterior temperature 
equal to the cruise exterior temperature. The equation is as follows: 


h^e/©H - Kt f (T s -T H )/L ss 2 rO 


For our design, the insulation must keep the hydrogen fuel at - 
4352F to ensure that it will remain in its liquid state. Furthermore, 
the cruise exterior temperature is approximated to be 120C£F. The 
material that was chosen for the insulation is Silica Fibers. In the 
calculation of the insulation thickness, the thickness of insulation 
was found to be 0.34 ft. The weight of this insulation is 
approximately 450 lbs. 

The thickness of insulation can be reduced further if a better 

/ 

insulation is used. One such insulation is Quartz Fiber which has a 
lower density (.66 at 70C£F), and a greater maximum temperature 
limit (25002F). The reason we have not included such fiber in the 
design was that we were unable to ublairi its thermal cunslanls. 



We strongly believe that there are better insulation materials 
other than those mentioned above. If further information on 
insulation materials can be obtained, the thickness of insulation will 
reduced substantially. 



Cost Analysis 
Thomas Greiner 


The cost analysis method used in determining the total cost of the 
hypersonic aircraft design was found in a book by Nicolai, Fundamentals of 
Aircraft Design (see references) and refered all dollar amounts to the cost 
in the year 1970. The final cost estimate was multiplied by the inflation 
factor of 2.35 to adjust for the current year 1990. The analysis based the 
cost on three main parameters, 

1) AMP.R Weight which is the weight of the dry aircraft minus the 
weight of the engines, the starter, all cooling fluids, wheels, brakes, 
instruments, auxiliary power and batteries. The value used was an average 
of the gold design team and the WAATS program weight estimate, 
A M ,PR= 18526 lbs. 


2) S which is the maximum speed at the "best" altitude in knots. 
This description is vague therefore the velocity for Mach 6 at 80,000 feet 
was used due to the large available thrust, 5=5178.22 knots 

3) Q which is the combined number of test aircraft and the number 
to be built during production. Since there is to be only one test aircraft to 
be built 0=1. 


Using the equation supplied in the book the cost of engineering hours 
was the largest cost at 1 7 billion dollars, see figure cost! while the 
smallest expense is the manufacturing and material cost at 1 1 million 
dollars. The materials estimate appears to be too low due to the unique 
alloys needed to cool and maintain structural integrity in the intense heat 
of hypersonic flight. 

The total price of the aircraft is 3.3 billion dollars which is most 
likely a gross underestimate. The intial development costs of the B-2 



Stealth bomber was approximately 40 billion dollars. Although the 
hypersonic aircraft would not need the technology of the B-2 bomber it 
stands as a good indication of the expense of applying new technologies to 
aircraft design. This aircraft would use the state of the art technologies 
in structural composites, ceramics in cooling, and propulsion to acheive 
its mission requirements. This total cost would be spread over a period of 
approximately five years. 



WIND TUNNEL MODEL 
by 

Bob stonebraker 

A 1/65 scale model of the preliminary design was built for 
tunnel testing and to aid presentations. This scale was selected 
to match that of other designs for relative comparison and to avoid 
overloading the tunnel measuring device with too large of a wing 
area. The model was constructed of wood and finished with black 
pigmented lacquer. Pine was chosen over balsa or fiber-glassed 
styrofoam for reasons of better durability and lower cost. Lacquer 
was chosen as the finish to allow polishing and provide a smooth, 
low drag surface. Construction took place only after approval of 
the preliminary three-view drawing. 

The model will be used to study the longitudinal stability and 
landing characteristics of the design. Tunnel testing will yield 
lift, drag and pitching moment for various angles of attack. 
Reynolds numbers up to 1.8 million can be obtained in the available 
tunnel with reference to the longitudinal dimension of the vehicle. 
Since this is a preliminary design, it is subject to change. Data 
from this model will be used in an additional design iteration and 
another model will be built and tested to incorporate changes. 



Conclusion 


This report has described and presented many concepts on how and why 
a hypersonic aircraft should be designed similar to the gold design for the 
given mission requirements. 

This paper attempted to find answers to questions concerning the 
areas of aerodynamics, thermodynamics, structural design, and aircraft 
configuration in a hypersonic enviroment. It can be concluded that these 
major design concerns have been addressed but the aircraft still remains 
in the early stage of development. 

Some problems as inlet design, fuel reduction, and dynamic stability 
still need to be resolved for the hypersonic aircraft. These problems 
currently exist within the aerospace industry and will have to be solved 
before hypersonic flight to be commonplace 
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Ip&flLE ! T /V / f LJ- &HJ P P Fa 

D Based On M=6 at 75,000 ft. 
a- 967.9481 v- 5B07.658B078 □- 1B5B.795 

CLinB ANGLE-4 DEGREES 

Uint 1404. B43 


fl- . 8-774 . 353FT/S 


Altitude 

CFt.) 

Temp . 
CR) 

Density 
C slug/Ft *3) 

Q U 

CSLUG/S*8 CFT/S) 

M 

40000 

389.93 

0.0005B7 


774.3537 

0.8 

40000 

3B9.99 

0 . 0005B7 


967.9481 

1 

40000 

389.99 

0.0005B7 


1161.530 

1 .8 

40000 

3B9.99 . 

0.0005B7 


1451.913 

1.5 

40000 

389.99 

0.000587 


1935. BB4 

8 

40000 

3B9.93 

0.0005B7 


8419.855 

2.5 

40000 

389.99 

0.0005B7 

1858.75 

8515.976 

8.599304 

45000 

383.99 ' 

0.000468 

1858.75 

8B35.813 

8.929734 

50000 

389.99 

0.000364 

1858.75 

3196.159 

3.302015 

55000 

389.99 

0 . 0002B7 

1B5B.75 

3608 .035 

3.781333 

60000 

389.99 

0.000886 

1B5B.75 

4059.854 

4.193695 

65000 

389.99 

.0.000178 

1858.75 

4574. 234 A 

4.785731 

70000 

389.99 

0.000140 

1B5B.75 

5154.304 

5.385013 

75000 

389.99 

0.000110 


5B07.651 

5.999999 

80000 

389.99 

0 .'000087 


5807.651 

5.999999 

B5000 

394.32 

0.00006B 


5839. B03 

5.999999 

30000 

408.48 

0.000053 


5B99.91B 

5.999999 

95000 

410.64 

0.000041 


5959.426 

5.999999 

100000 

418.79 

0.000038 


6018.874 

5.999999 

100000 

418.79 

0.000038 


6519.798 

6.5 

100000 

418.79 

0.000038 


7081.381 

7 

100000 

41B.79 

0.000038 


7582.844 

7.5 

100000 

410.79 

0.000038 


8084.367 

0 

100000 

41B.79 

0.000038 


B585.B90 

8.5 

100000 

41B.79 

0.000038 


9087.413 

9 

100000 

418.79 

0.000038 


9528.936 

9.5 

100000 

418.79 

0.000038 


10030.45 

10 

100000 

418.79 

0.000038 


10030.45 

10 


5s 




/^lT 


Mach # 

^Altitude 

Maos <flew 



Y 

' •• 

• (ft) 

(lb/sec) 

lIBlllifl 


(f t**2 

.80 

40000.0 

495.0 

33.804 * 

33. 804 

4.829 

1.00 

40000.0 

495.0 

27.043 

27.449 

3.921 

1.20 

40000.0 

495.0 

22.536 

23.077 

3.297 

1.50 

40000.0 

478.0 

17.410 

18. 002 _ 

2.572 

2.00 

40000.0 

380. 0 

10. 380 

10.8997. 

1.557 

2.50 

40000. 0 

304.0 

6. 643 

7. 108/ 
6. 563 

1.015 

2.60 

40000.0 

290.0 

6.094 

.938 

2.93 

45000.0 

250.0 

5.922 

6.467 

.924 

3. 30 

50000.0 

210.0 

5.611 

6. 200 A 

.886 

3.72 

55000.0 

175.0 

5.268 

5. 889 

.841 

4. 19 

60000.0 

148.0 

5.023 

5.681 

.812 

4.72 

65000. 0 

126.0 

4.821 

5.532 

.790 

5.32 

70000.0 

110.0 

4.741 

’ 5.529 

.790 

6.00 

75000 . 0 

94.0 

4 . 560 

5.416} 

.774 

6.00 

80000.0 

94.0 

5.789 

6.874. ' 

.982 

6 . 00 

85000. 0 

94.0 

7.379 

8.763 7 

1.252 

6.00 

90000. 0 

94.0 

9.419 

11. 185 

1.598 

6. 00 

95000. 0 

94.0 

1 1.956 

14. 198 

2. 028 

6. 00 

100000. 0 

94.0 

15. 106 

17.939 

2.563 
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Table IN3 : 

RAMP LENGTHS VS. 

MACH NUMBER. 

M=6 , NR=8 

M=5, NR=8 

M=4, NR=8 

7.2746 

5.6746 

3.6746 

3.9430 

3.4430 

2.0000 

2.0864 

2.1000 

1.6000 

1.0900 

1.1000 

1.8430 

M=3 , NR=5 

M=2, NR=3 

M= 1 .4, NR= 

3.6746 

1.6746 

- 

2.0000 

- 

- 



tMlb ; 


tl££f\ 

Mach # = 

. \t§ llp^' 

1.40 - 3 

• 

THETA 

BETA 

M ( I ) 

M < i+i ) 

POR , PR 

TR 

XT 

5.25 

53.23 

1.400 

1.205 

.99815 1.3005 

1 . 0785 

2.80 

TPOR 

.99815 

TPR 
1 . 3005 

TTR 

1 . 0785 

TPOTR TPTR 

. 99043 1 . 9881 

‘rvAnKftnn W R n n A 

TTTR 
1 . 2203 

A A A aaa A 

MNT 

.8388 


rtme: 

1SG 

THtTfi, 

Mach # 

B€ 7V t 
= 2.00 . 

B rriciE/vtyfo 

n 

THETA 

BETA 

M ( I ) M ( 1+1 ) 

POR 

PR 

TR 

XT 

5. 25 

si!' A" m Jl. 

'2.000 1.812 

. 99758 

1 M 0«n!>vL'« 1 

1 . 0864 

5.45 

5. 16 

38.32 

1.812 1.634 

. 99806 

1 . 3063 

1 . 0799 

4.75 

10.41 

50.28 

1.634 1.261 

. 98607 

1.6773 

1 . 1 639 

3.12 
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A 

****************************************** 

TPOR 

TPR 

TTR 

TPOTR 

TPTR 

TTTR 

MNT 

.98178 

2. 9209 

1 . 3655 

.96753 

4.9344 

1 . 5928 

. 8064 


~ fi-tfLF ; J A/ y 

7T/£ TA BBT A, 

Mach # = 3.00 

ctj f 0 r- 


THETA BETA 

M ( I ) 

M <1 + 1 ) 

POR 

. PR 

TR 

XT 

e* ncr ot tt 

iJ O ■ O O 

3. 000 

2.738 

.99388 

1 . 4803 

1. 1206 

8. 35 

10.13 29.60 

2.738 

2.282 

.96924 

1 . 9669 

1.2241 

6.34 

4.27 29.36 

2.282 

2.4 16 

. 99826 

1 . 2938 

1 . 0769 

6. 40 

5.90 33.16 

2. 1 16 

1 . 899 

.99621 

1 . 3966 

1. 1013 

5.51 

13.75 46.23 

1.899 

1 . 40 1 

. 96506 

2.0271 

1 . 2362 

3. 45 




*********#*******•****************•*•*••****■**.* 

TPOR TPR 

TTR 


TPOTR 

TPTR 

TTTR 

MNT 

.92452 10.6647 

2.01 1 1 


.88575 

22. 6281 

2. 5240 

. 7395 
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AfitE ; 

r/v 7 

THE7A, TEFF ICT£ wcy jor / 

Mach # = 4.00 

THETA 

BETA 

M < I ) 

M(I'+1 > 

FOR 

PR 

TR XT 

5- 25 

18. 22 

4 . 000 

3.621 

.98701 

1 . 6573 

1.1 596 1 0 . 64 

4.41 

19. 15 

3. 621 

3. 342 

. 99393 

1 . 4787 

1.1202 10.08 

3. 43 

19.82 

3. 342 

3. 145 

,/99763 

1 . 3306 

1 . 0858 9.71 

5. 42 

22.50 

3. 145 

2.859 

. 99244 

1 . 5236 

1.1303 8.. 45 

1 1 . 04 

29.42 

2.859 

2 . 340 

. 957 1 0 

2. 1343 

1 nr -/ cr a \ 

X a v_j / W • -L. J. 

3 . 55 

28.04 

2.340 

2. 199 

. 99893 

1 , 2452 

1.0650 6.57 

12. 25 

38. 15 

. 2. 199 

1 . 735 

. 96790 

1 . 9866 

1.2281 4.46 

13. 75 

51.72 

i a / o 5 

1.229 

. 96715 

1 . 9974 

1.22.02 . 2.76 

TPOR 

TPR 

TTR 


TPOTR TPTR TTTR MNT 

.86929 

52. 3936 

3.2256 


. 86032 

83.6254 

3.6975 -.8245 



77/5 7 ' A y /IE 7 A y EFFJLcEc. 

A /C ^ r ’• 

~mle ; 



Mach # 

= 5 . 00 



THE T A 

BETA 

M < I ) 

M ( I + 1 ) 

FOR 

PR 

TR XT 

5 . 25 

1 5 • 27 

5 . 000 

4.469 

. 97634 

1 . 8565 

1.2015 12.82 

6. 08 

1 7 . 32 

4.469 

3.954 

. 97370 

1 . 8989 

* 1 .2103 11.22 

7 . 16 

19.94 

3. 954 

3. 448 

. 97013 

1 . 9538 

1.2214 9.65 

9 . 0 1 

23 - / / 

3. 448 

2.921 

. 96073 

2.0865 

1.2481 7.95 

8. 13 

26.30 

2.921 

2.529 

. 98034 

1 . 7879 

1 . 1 873 7 . 08 

9 „ 63 

3 i - i 8 

2.529 

2. 127 

. 977'68 

1 . 8341 

1.1969 5.78 

i 2 « 25 

39.39 

2. 127 

1 . 671 

. 96985 

1 . 9579 

1.2223 4.26 

X » / wJ 

54.75 

1 . 67 1 

la i O X 

. 96663 

2.0049 

1.2318 2.47 



— 
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* * -S- * * -S- * ¥r * ¥; * -*■ ¥r *■ * ###•*•*■*••*•*"**•*••*• -S • -R * # Sr 

TPOR 

TPR 

TTR 


TPOTR 

TPTR 

TTTR MNT 

- 7 9 6 1 5 

184. 9801 
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. 79345 
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ZPf 74, 
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for; 


7 /v 7 

( 

Mach # 

~ 6.00 



THETA 

BETA 

M (I ) 

M ( I + 1 ) 

POP 

PR 

TR XT 

5,. 25 

13. 36 

6. 000 

5.282 

„ 961 46 

2.0766 

1.2461 14.73 

6. 25 

15,. 47 

5.282 

4 . 594 

. 95595 

2. 1491 

1.2604 12.65 

t •*:« c: 

17.94 

4.594 

3 . 956 

. 95443 

2. 1686 

1.2643 10.81 

8 . 50 

21 . 07 

3 956 

3. 357 

. 95252 

2. 1926 

1.2690 9.08 

9 * 7 5 

24 „ 91 

3 . 357 

2. 806 

. 95462 

2. 1661 

1 . 2638 7 . 54 

11. 00 

29 . 82 

2 . 806 
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. 95932 

2. 1052 

1 . 2518 6. 1 1 
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56. 6,5 

2.299 

1 . 822 

96497 

2. 0283 

1.2365 4.71 

i 3 - 75 
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1 . 822 . 

1 . 323 

.96644 

2. 0076 
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TABLE ADI Components of the drag coefficient for the wing at different 
Mach numbers. 


M 

C Df 

C Do 

c Dle 

^[)w 

0.8 

0.005156 

- 

- 

- 

0.9 

0.005156 

0.001643 

- 

- 

1.0 

0.005156 

0.013970 

- 

- 

1.1 

0.005026 

0.014010 

- 

- 

1.2 

0.004872 

0.014092 

- 

- 

1.5 

0.004382 

- 

0.000612 

0.01263 

2.0 

0.004021 

- 

0,001225 

0.00545 

2,5 

0,003531 

- 

0.001904 

0,00510 

3.0 

0.003196 

- 

0.002522 

0,00508 

3,5 

0.002835 

- 

0.003021 

0.00520 

40 

0,002578 

- 

0.003398 

0.00529 

4.5 

0.002320 

- 

0.003675 

0.00534 

5.0 

0.002114 

- 

0,003877 

0.00537 

5.5 

0,001907 

- 

0,004066 

0.00542 

6,0 

0 001730 

- 

0 004135 

0,00537 

6.5 

0,001624 

- 

0,004218 

0.00536 

7.0 

0,001443 

- 

0 004280 

0,00534 

7,5 

0.001392 

- 

0,004329 

0,00531 

8.0 

0,001289 

- 

0 004366 

0,00529 

8,5 

0,001 160 

- 

0,004400 

a Aftrm 

u.Ovoz / 

9.0 

0,001 108 

- 

0 .004420 

0.00524 

9,5 

0.001057 

- 

0,004440 

0.00521 

1 0 .0 

0,000830 

- 

0.004540 

0.00519 



TABLE AD2 Components of the drag coefficient for the body at different 
Mach numbers. 


n 

C Df 

C Db 

C Dp 

C Dp 

C DN 2 

C DA 

C DA ( NC ) 

0.8 

0.08808 

0.000587 

_ 

- 

- 

- 

- 

0,9 

0,08007 

0,006136 

0,00795 

- 

- 

- 

- 

1.0 

0.08007 

0,016674 

0,00795 

0.063 

- 

- 

- 

1,1 

0,07807 

0,012272 

0.00375 

0145 

- 

- 

- 

1,2 

0,07567 

0.009972 

0.00100 

0,154 

- 

- 

- 

1,5 

0. 06806 

0,001022 

- 

- 

0,10869 

0,03168 

0 0222 

2.0 

0,06245 

0.002315 

- 

- 

0,10640 

0,02620 

0,0222 

2,5 

0 05485 

0,005676 

- 

- 

0,09496 

0,02200 

0.0222 

3.0 

0.04964 

0,007670 

- 

- 

0.08695 

0.02020 

0,0222 

3.5 

0.04404 

0,007593 

- 

- 

0.08009 

0,01843 

A A A AAA 
V i 

4,0 

0.04004 

0,007287 

- 

- 

0.07780 

0,01613 

0.0222 

4,5 

0.03603 

0,006443 

- 

- 

0,07440 

0,01469 

0,0222 

5,0 

0.03283 

0,005753 

- 

- 

0.07322 

0.01354 

0-0222 

5,5 

0. 02963 

0,004909 

- 

- 

0,07151 

0.01267 

0.0222 

6,0 

0,02722 

0.04295 

- 

- 

0.06980 

0,01 152 

0.0222 

6.5 

0.02522 

0.00399 

- 

- 

0.06868 

0.01037 

0 . 0222 J 

7,0 

0.02242 

0.03990 

- 

- 

0.06750 

0.01008 

0.0222 

7,5 

0.02162 

0.03990 

- 

- 

0.06693 

0.00922 

0.0222 

8,0 

0,02002 

0.03990 

- 

- 

0.06636 

0.00864 

0.0222 

8.5 

0.01842 

0.03990 

- 

- 

0.06579 

0.00306 

0.0222 

9.0 

0.01722 

0.03990 

- 

- 

0.06521 

0.00777 

0,0222 

9.5 

0.01644 

0.03990 

- 

- 

0.06464 

0.07488 

0.0222 

10.0 

0.01353 

0.03990 

- 

- 

0.06407 

0.006912 

0.0222 



TABLE AD3 Cq 0 for the different Mach numbers. 


li 

^ c 0oHotal 

0.0 

0.012051 

0.5 

0.012938 

0.8 

0,012988 

0.9 

0.014719 

1.0 

0,034345 

1.1 

0.039220 

1.2 

0.039493 

1.5 

0.037375 

2.0 

0.028804 

2.5 

0.025900 

3.0 

0.024200 

3.5 

0.022705 

4.0 

0.021767 

4.5 

0.020770 

5.0 

0.020104 

5.5 

0.019427 

6.0 

0.018758 

6.5 

0.018237 

7.0 

0.017658 

7.5 

0.017373 

8 0 

0.017012 

8.5 

0,016643 

9,0 

0.016390 

9,5 

0.016125 

10,0 

0,015517 





TABLE AD4 Shows Cj_ for the different Mach numbers, 


M 


1.5 

- 

0.4292 

- 

5.7 

2.188 

2.0 

- 

0.6650 

- 

5.45 

2.092 

2.5 

- 

0.8794 

- 

5.0 

1.9193 

3.0 

0.921 

- 

4.70 

- 

1.662 

3.5 

0.777 

- 

4.60 

- 

1.371 

4.0 

0.673 

- 

4.40 

- 

1.136 

4.5 

0.594 

- 

4.30 

- 

0.980 

5.0 

0.532 

- 

4.22 

- 

0.861 

5,5 

0.482 

- 

4.18 

- 

0.773 

6.0 

0.440 

- 

4.16 

- 

0,703 

6.5 

0.406 

- 

4.14 

- 

0.645 

7.0 

0.376 

- 

4.14 

- 

0.598 

7,5 

0 350 

- 

4.13 

- 

0,556 

0.0 

0.328 

- 

4.13 

- 

0.520 

0.5 

0.309 

- 

4.12 

- 

0 488 

9.0 

0.291 

- 

4.10 

- 

0,458 

9.5 

0.276 

- 

4.07 

- 

0.431 

10,0 

0.262 

- 

4.05 

- 

0.407 



TABLE AD5 


M 


0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

1.5 
2.0 

2.5 

3.0 

3.5 

4.0 

4.5 

5.0 

5.5 

6.0 

6.5 

7.0 

7.5 

8.0 

8.5 
9.0 

9.5 

1 0 0 


L /rad, Cj_ /drqree, and K for different Mach numb 


C L /rad 

K 

1.633 

0.31032 

1 .640 

0.31032 

1.645 

0.31032 

1.651 

0.31032 

1.656 

0.31032 

1.662 

0.31032 

1.670 

0.31032 

1.689 

0.31032 

2.188 

0,35000 

2.092 

0.43000 

1.919 

0.57100 

1.662 

0.60200 

1.371 

0.72940 

1.136 

0.88020 

0.980 

1 .02000 

0.861 

1.16100 

0.773 

1 .29400 

0.703 

1 .42200 

0.645 

1 .55000 

0.598 

1 .67200 

0.556 

1.79900 

0.520 

1 .92200 

0,488 

2.04900 

0.458 

2.18300 

0.431 

2.32000 

0.407 

2.45700 



TABLE AD6 

Cl, Cq, and C^ 

/C D for 

our mission profile. 

M 

C L 

C D 

c L / c D 

0.8 

0.37966 

0.05902 

6.43292 

1 

0.24244 

0.05541 

4.37519 

1.2 

0.16807 

0.05135 

3.27296 

1.5 

0.10727 

0.05075 

2.11384 

2.0 

0.06003 

0.03755 

1.59855 

3.0 

0.. 03597 

0.03103 

1.15919 

4.0 

0.03542 

0.027319 

1.2964 

5.0 

0.03524 

0.02657 

1.32637 

6.0 

0.03531 

0.02522 

1.40017 

7.0 

0.03371 

0.02397 

1.40635 

8.0 

0.3308 

0.02337 

1.41544 

9.0 

0.03108 

0.0226 

1.35729 

10.0 

0.02883 

0.02144 

1 .34487 

10,0 

0.02947 

0.02153 

1 .36874 

10.0 

0.03010 

0.02162 

1 ,39206 

10.0 

0,03074 

0.02172 

1,4154 

10.0 

0,03139 

0,02182 

1,4388 



MACH .B 

Cm eg 

AnglB of Attack 

CL 

CD 

. .0009005 

0.0 

0.0 

.0129005 

- . 04B9BB 

3.459 

0.1 

.016090 

-.094403 

0.917 

0.2 

.025400 

-.141400 

10.370 

0.3 

.040920 

-.107994 

13.035 

0.4 

.002640 

-.234109 

17.293 

0.5 

.090559 

-.279044 

20.752 

0.0 

.124705 

-.325120 

24.211 

0.7 

.165047 

-.309902 

27.509 

0.9 

.211596 

-.414405 

31.129 ‘ 

0.9 

.264351 

-.450390 

34.500 

1.0 

.323312 

Table 51 : 

Vs. ai, C L , ih=,t) 




MACH 2 



Cm eg 

Angle of Attack 

CL 

CD 

.0019975 

0.0 

0.0 

.020004 

-.0457407 

2.553 

0.1 

.033574 

-.0920155 

5.125 

0.2 

, 04B2B4 

-.1392000 

7.509 

0.3 

.072634 

-.1049205 

10.252 

0.4 

.106724 

-.2299090 

12.015 

0.5 

.150554 

-.2743300 

15.370 

0.6 

.204124 

-.3100274 

17.941 

0.7 

.267434 

-.3010434 

20.504 

0.B 

.340404 

-.4033039 

23.057 

0.9 

.423274 

-.4450409 

25.630 

1.0 

.515004 

Table 52 ' 

v?, at, £<.,<„ in-t) 




MACH 10 



Cm eg 

Angle of Attack 

CL 

CD 

.0010701 

0.0 

0.0 

.015517 

-.003715 

1.400 

0.01 

.01576 

-.000472 

2.015 

0.02 

.01650 

-.013195 

4.223 

0.03 

.01773 

-.017004 

5.631 

0.04 

.01945 

-.022539 

7.039 

0.05 

.02166 

-.027150 

0.447 

0.06 

.02436 

-.031745 

9.054 

0.07 

.02756 

-.030299 

11.262 

0.00 

.03124 

-.040017 

12.670 

0.09 

.03542 

-.045302 

14.070 

0.10 

.04009 

-.000272 

20.155 

0.20 

.11390 

-.127034 

42.233 

0.30 

.23665 

-.103909 

56.310 

0.40 

. 40B64 

-.195735 

70.300 

0.50 

.62977 

-.225074 

04.460 

0.50 

.90004 


Table S3: (- WCi) Vs. c \ , U , C 0 (M-IO) 
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Uehicle Summary 


Uehicle Name: tscram 


:le Discription 

: UHANNED 

SCRAI'IJET TESTBED 



File = B : GOLD . 

DAT 

Output 

File 

= B: GOLD. OUT 

Date = 01-01- 

1380 


Date 

= 01-01-1380 

Time “ 00:01: 

27 


Time 

- 00:03:13 

h •= 73 

Height " 

8 Span 

32 

TRoot =* 5.1 

y = 1,570 

SWing = 

508 SUert = 

51 

SHorz = 0 

er of : Ramjets 

- 5 

Rockets = 0 

TurbRamJets = 0 

st of: - 

100,000 

0 


0 

10,000 


Fuel Density = 4,383 


LOX - 

ts: GT0W = 

44,133 

Pauload = 

0 

Dry = 33,S03 

Landing = 

34,013 

Entry = 34,033 

s: GTOW-'S - 

86 . 34 

AR = 1.380 T/GTOW = 

11.33 

T/S = 385.00 


Weight Statement 

Group i: ^ero surfaces 

Wing 
Uert ical 
Horizontal 
Fairing 

Group 2: Eody structure 

Easic body 

Secondary 

Thrust 

Integral Fuel tanks 
Integral Ox tanks 


Group 3: Thermal Protection System 

Uehicle insulation 
Cover panels 


Group 4: Launch and Recovery Gear 

Launch gear 
Landing gear 


Group 5: Propulsion 

Rocket engines 
Ram J ets 
Turboramjet 

Nonstructura 1 Fuel tank 
Nonstructural Ox tank 
Fuel tank insulation 
Ox tank insulatiqn 
Fuel system 
Oxidizer system 
Pressurization system 
Inlets 


310 

543 

361 

0 

0 


6, 1B2 
1,533 
3, 134 
1,306 
0 


2,352 

0 


110 

1,527 


0 

7,000 

0 

0 

0 

704 

706 

650 

233 

1,025 

175 


12, 821 

2,352 
1 , 63B 
10,353 


FfeufiE WEIGHT 1 


13 


Weight statement Far:tscram Page 2 


Group 

G: 

Orientation Control System 


B34 



Engine gimbal system 

0 




Attitude control system 

1B3 




Aerodynamic controls 

573 




Seperation system 

132 




ACS tankage 

0 


Group 

8: 

Power supply 


754 



Electrical System 

716 




Hydraulic/Pneumatic Sys 

36 


Group 

10: 

Avionics 


3,153 

Group 

14: 

Crew Provisions 


422 



Uehicle Dry 

Weight ' 

33,603 

Group 

17: 

Crew 


35 

Group 

18: 

Payload 


0 

Group 

21 : 

Residual Propellant 


75 



Trapped fuel 

75 




Trapped Oxidizer 

0 




Landing 

Weight 

34,013 

Group 

22: 

Reserve Propellants 


BO 



Fuel 

BO 




Oxidizer 

0 




ACS fuel 

0 




ACS oxidizer 

0 




Entry 1 

height 

34,033 

Group 

23: 

Inflight Lasses 


40 



Fuel 

40 




Oxidizer 

0 


Group 

25: 

Main Propellants 


10,000 



Fuel 

10,000 




Oxidizer 

0 




Qgoss 

Weight 

44,133 


FI6.U*E WEIOtri 

c.owtiwue.0 



Gold Weight Estimate 
Thomas Greiner 


Structure- 


Skin 

7955 lb. 

Internal 

2386 lb. 

Total ' 

10,341 lb. 

Engines 


Scram jets 

1500 1b. x 4 12,000 lb. 

T C Q 1 

6500 lb. 6500 !b. 


Total 18,500 1b. 

Electrical 


Computers 

75 lb. 

Batteries 

300 lb. 

Avionics 

100 lb. 

Transmitter 

55 lb 

Total 

510 lb 

Cooling Systems 


Cedant 

315 lb 

i r *• r* \ i • t y ■*' r. 
i • J •. < 'UU*,lVi‘ 

435 lb. 

U Ci t 

• . w Av, i fOi 

• 485 1b. 

Pumps 

75 lb. 

Total 

1310 lb 

Fuel Tam- 

4013 lb 

insulaticn 

400 lb. 

Landing Gear 

200 lb. Empty Weight 36 : 


Fuel i 1 ,500 lb. Gross Weight 47,774 lb. 

Figure: Weight2 ORiGSNAL PASS m 

OF POOR QUALITY 




Rocket Booster Fuel Estimate for Hypersonic Vehicle 


Aircraft weight 47000 lb 
Estimated booster weight 15000 lb 
Estimated total weight 62000 lb. 



Specif ic 

Impulse 240 

s 







Ig acceleration 





Mach 

Vel 

(ft/s) 

Thrust 

(lb) 

t 

(S) 

Total 

t 

(s) 

W-dot 

(lb/s) 

Rocket 

Fuel 

(lb) 

SFC 

Seconds 

0.8 

774.35 

0 

0 

0 

0 


0 0.000099 


1 967.94 62000 6.012111 6.012111 258.3333 1553.128 0.0002 
1.2 1161.53 60436.18 6.012111 12.02422 251.8174 1513.954 0.000196 
1.5 1451.913 58888.26 9.018105 21.04232 245.3677 2212.752 0.000194 


2 1935.884 56610.25 15.03015 36.07248 235.8760 3545.254 0.00021 
2.5 2419.85 52887.42 15.03 51.10248 220.3642 3312.075 0.000219 

12137. 16 


Figure: Booster2 



TFRJ 

TFRJ 

A/C 

rocket 

Section 


Fuel 

Fuel 

Weight 

Volume 

length 

Drag 

<lb/s) 

(lbs. ) 

(lb- ) 

(ft~3) 

ft 


7177 0.710523 0 62000 000 

8886.99 1.777398 10.68591 60436.18 15.30176 3.117252 74.55013 

14413.38 2.825022 33.96870 58888.26 14.91580 3.038626 71.21642 

15984 3.100896 65.25007 56610.25 21.80051 4.441168 103.0257 

23442 4.92282 177.5783 52887.42 34.92861 7.115604 178.6808 

31258 6.845502 349.8221 49225.52 32.63128 6.647595 174.0826 


637.3051 24.36024 601.5558 


Figure: Booster2 continued. 


IS 



Q based on M = 6 at 70,000 


a= 967.9421 v= 5807.6528078 

CLIMB ANGLE=4 DEGREES 

M=. 8=774. 353FT/S 
Altitude Temp. Density 

(ft.) ( R) (slug/ft~3) 

Q= 

Q 

( SLUG/S' - 

2359.837 

Vint 

V 

2 (FT/S) 

1404.843 

M 

40000 

389.99 

0.000587 


774.3537 

0.8 

40000 

389.99 

0.000587 


967.9421 

1 

40000 

389.99 

0.000587 


1161.530 

1.2 

40000 

389.99 

0.000587 


1451.913 

1.5 

40000 

389.99 

0.000587 


1935.884 

2 

40000 

389.99 

0.000587 


2419.855 

2.5 

40000 

389.99 

0.000587 

2359.837 

2834.895 

2.928786 

45000 

389.99 

0.000462 

2359.837 

3195.274 

3.301100 

50000 

389.99 

0.000364 

2359.837 

3601.297 

3.720571 

55000 

389.99 

0.000287 

2359.837 

4058 . 621 

4.193041 

60000 

389.99 

0.000226 

2359.837 

4573.796 

4.725278 

65000 

389.99 

0.000178 

2359.837 

5154.053 

5.324754 

70000 

389.99 

0.000140 

2359.837 

.5807.651 

5.999999 

75000 

389.99 

0.000110 


5807.651 

5.999999 

80000 

389.99 

0.000087 


5807.651 

5.999999 

85000 

394.32 

0.000068 


5839.803 

5.999999 

90000 

402.48 

0.000053 


5899.918 

5.999999 

95000 

410.64 

0.000041 


5959.426 

5.999999 

100000 

418.79 

0.000032 


6018.274 

5.999999 

100000 

418.79 

0.000032 


6519.798 

6 . 5 

100000 

418.79 

0.000032 


7021.321 

7 

100000 

418.79 

0 . 000032 


7522.844 

7.5 

1000Q0 

418.79 

0.000032 


8024 . 367 

8 

100000 

418.79 

0.000032 


8525.890 

8.5 

100000 

418.79 

0.000032 


9027.413 

9 

100000 

418.79 

0.000032 


9528.936 

9 . 5 

100000 

418.79 

0.000032 


10030.45 

10 

100000 

418.79 

0.000032 


10030.45 

10 


F igure : Prof i le 1 


■?q 



Total Total Total 

dist t t x Dist Dist acc G 

(angle 4 )( seconds )( seconds ) (ft) (ft) (miles) ft/s'"2 



5 

5 

3871.768 

3871.768 

0.733289 


0 

0 


20 

25 

17422.95 

21294.72 

4.033092 

9 . 

679421 

0.300603 


20 

45 

21294.72 

42589.45 

8.066184 

9 . 

679421 

0.300603 


20 

65 

26134.43 

68723.89 

13.01588 

14 

.51913 

0.450904 


20 

85 

33877.97 

102601.8 

19.43217 

24 

. 19855 

0.751507 


20 

105 

43557.39 

146159.2 

27.68167 

24 

. 19855 

0.751507 

71677.99 

27.28121 

132.2812 

71505.44 

217664.7 

41.22437 

15 

.21341 

0.472466 

71677.99 

23.77312 

156.0543 

71505.44 

239170.1 

54.76707 

15 

. 15907 

0.470778 

71677.99 

21.09239 

177. 1467 

71505.44 

360675.6 

68.30977 

19 

.24974 

0.597818 

71677 . 99 

18.71507 

195.8618 

71505.44 

432131.0 

81.85247 

24 

.43610 

0.758885 

71677 . 99 

16.60670 

212.4685 

71505.44 

503686.4 

95.39516 

31 

. 02211 

0.963419 

71677.99 

14 . 73665 

227.2051 

71505.44 

575191.9 

108 . 9378 

39 

. 37510 

1.222829 

71677.99 

13.07789 

240.2830 

71505.44 

646697 . 3 

122.4805 

49 

. 97733 

1.552091 

71677.99 

12.34199 

252.6250 

71505.44 

718202.8 

136.0232 


0 

0 

71677.99 

12.34199 

264 . 9670 

71505.44 

789708.2 

149.5659 


0 

0 

71677 . 99 

12.30792 

277.2749 

71505.44 

861213.7 

163. 1086 

2 . 

612280 

0.081126 

71677.99 

12.21119 

289.4861 

71505.44 

932719.1 

176.6513 

4 . 

922911 

0. 152885 

71677.99 

12.08301 

301.5741 

71505.44 

1004224 . 

190.1940 

4 . 

922911 

0.152885 

71677 . 99 

11.96857 

313.5427 

71505.44 

1075730 

203.7367 

4 . 

916878 

0. 152898 


43 

356.5427 

269568.5 

1345298 . 

254 . 7914 

11 

.66334 

0.362215 


43 

399.5427 

291134.0 

1636432 . 

309.9304 

11 

.66332 

0 . 362215 


70 

469.5427 

509045.7 

2145478 . 

406 . 3406 

7 . 

164613 

0 . 222503 


70 

539.5427 

544152.4 

2689630 . 

509 . 399 ? 

7 . 

164613 

0.222503 


70 

609.5427 

579259.0 

3268689 . 

619. 1079 

7 . 

164613 

0.222503 


70 

679.5427 

8 14365 . 6 

3883255 . 

735.4650 

7 . 

164613 

0.222503 


70 

749.5427 

649472.2 

4532727 . 

858.4711 

7 . 

164613 

0.222503 


70 

819.5427 

634578.3 

5217306 . 

988 . 1262 

7 . 

164613 

0.222503 


120 

939 . 5427 

1203655 

6420961 . 

1216.091 


0 

0 


M=.8 to 6 

time minutes x-dist miles 


939 . 542 ? 15.65904 5217306 . 983.1262 


M=6 to 10 

time minutes x-dist miles 


626 10.43333 5345231 . 1012.354 


F igure : Prof iiel 



Thrust 

acc 

Drag 

Thrust 

req 

SFC 

hour 

0 

5977 

5977 

0.357 

13526.25 

7252 

20778.25 

0.721 

13501.23 

11630 

25131.23 

0.709 

20207.21 

12849 

33056.21 

0.7 

33582.08 

18795 

52377.08 

0.8 

33407.13 

25046 

58453.13 

0.78 

20883.09 

12849.4 

36871.53 

0.63 

20725.63 

11500.0 

35296.59 

0.7 

26220.89 

11000.0 

40280.48 

0.73 

33154.72 

10153.7 

46355.98 

0.75 

41916.43 

9000.0 

53951.39 

0.845 

52945.65 

10000.0 

65965.94 

1 

66782.81 

10158.0 

79942.26 

1.063 

0 

11300.0 

14279.91 

1 

0 

11300.0 

14276.50 

1.02 

3457.624 

10990.0 

17420.64 

1.03 

6506.606 

10750.0 

20225.34 

1.045 

6495.645 

10670.0 

20129.38 

1.045 

6476.899 

10652.0 

20087.71 

1.12 

15336.78 

10652 

25988.78 

2 

14878.42 

11233 

26111.42 

2 

8856.467 

11790 

20646.46 

2 

3499.162 

12292 

20791 . 16 

2 

3139.354 

12832 

20971.35 

2 

7776.427 

13406 

21182.42 

9 

7409 . 848 

14007 

21416 . 84 

2 

7038 . 647 

14708 

21746.64 

2 

0 

14703 

14703 

9 


fuel Volume 


Vn cj I a> 1 

SFC Ib/s Weight 

seconds burned 


0.000099 0.592719 2.963595 
0.000200 4.161421 83.22843 
0.000196 4.949456 98.98913 
0.000194 6.427597 128.5519 
0.000222 11.63935 232.7870 
0.000216 12.66484 253.2969 
0.000175 6.452518 176.0325 
0.000194 6.863227 163.1603 
0.000202 8.167986 172.2824 
0.000208 9.657497 180.7408 
0.000234 12.66359 210.3005 
0.000277 13.32387 270.0326 
0.000295 23.60517 308.7058 
0.000277 3.966644 43.95629 
0.000283 4.045009 49.92347 
0.000286 4.984240 61.34564 
0.000290 5.870969 71.69152 
0.000290 5.843113 70.63167 
0.000311 6.249511 74.79773 
0.000555 14.43821 1241.686 
0.000555 14.50634 1247.545 
0.000555 11.47025 1605.836 
0.000555 11.55064 1617.090 
0.000555 11.65075 1631.105 
0.000555 11.76801 1647.522 
0.000555 11.89824 1665.754 
0.000555 12.08147 1691.405 
0.000555 8.171111 1961.066 


A/C 

Weight 


45000 

44997.03 
44913.60 
44814.81 
44686.26 
44453.47 
44200. 18 
44024 . 15 
43860.99 

43688.70 

43507.96 
43297.66 

43027.63 
42718.92 

42669.97 

42620.04 

42558.70 
42487.01 
42416.37 
42341.58 
41099.89 
39852.34 
38246.51 
36629.42 
34998.31 
33350 . 79 
31685.03 

29993.63 


15006.36 3387.441 


F igure : Prof ilel 
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CDo VS. MACH 


FIG AD3 




CD vs. MACH 


FIG AD4 



MAC II 



CL vs. MACH 


FIG AD5 
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L/D vs. MACH 


FIG AD6 



MACII - 



CL alpha vs. MACH * 


FIG AD7 
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Cm c.g. vs. Alpha (Mach .8) 
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LEADING EDGE 
RADIUS REQUIREMENT 


COOLANT RETURN VIA 
UPPER PANEL 

COOLED END CAP 


COOLANT 
SPRAY BAR 


COOLANT RETURN VIA 
LOWER PANEL 



TITANIUM CONCEPT 


SUPER PLASTICALLY 
FORMED SKIN 


COOLANT 

PASSAGE 



DIFFUSION BONDED 
STIFFENER 


VIEW A*A 




TYPICAL COOLED PANEL CONFIGURATIONS 

1NGER-STIFFENEDSK1N HONEYCOMB SANDWICH INSULATED/COOLED PANEL 


STRINGER-STIFFENED 
PLATE-FIN SANDWICH 



POTENTIAL POTENTIAL 

AD VANTA GES DISADVANTA GES/L/M/TA TIONS 

• USE OF ALUMINUM MATERIALS • REQUIRES LARGE FUEL HEAT SINK 

• LOW WEIGHT • COMPLEX SYSTEM ARRANGEMENT 

• VOLUMETRIC EFFICIENCY • MAY REQUIRE ADDITIONAL THERMAL 

PROTECTION DEVELOPMENT 




V\C^OC«~ xV ' 1 



///// //// //// ///////i 


HEAT SHIEL 







n 2 purge 


, JHIGH TEMP INSULATION W& 


CLOSED CELL INSULATION^ 





Total Development and Evaluation 


Airframe Engineering 
Development Support 
Flight Test Aircraft 
Engines & Avionic 
Manufacturing Labor 
Material & Equipment 
Tooling 

Quality Control 
Flight Test Operations 
Test Facilities 


1,697,859,405.60 
1,384, 147,764.00 
501,706,296.40 

20.414.855.50 

106.997.571.50 
10,931,449.70 

350,190,800.20 

13,171,619.50 

32,259,410.80 

OOO 

Total $3,615,972,876.80 


Figure: Costl 





APPENDIX A 


WING; 

SUBSONIC: 

C Do =C Df = c f U + 2*t/c+100*(t/c) 4 ] S wet /S ref 
C Df : Skin friction drag coefficientwhich is constant in the 
subsonic region. 

C f : Turbulent flat plate skin friction coefficient, 

(Nicolai, Figure E2). 

TRANSONIC: 

^Do =l "Df + ^Do =< “f H + t/c] S we j-/S re f + ACp 0 

AC Do : Wave drag coefficient (Nicolai, Figure 1 1.10). 

SUPERSONIC: 

For wings with round-nosed airfoil sections and supersonic leading 
edge use : 

c Do“ c Df tC Dle*l 1 6/3/J3]Mt/c> 2 *S wet /S ref 

* s wet/S re f 
c r [C fc /C r 1 ]*C fJ 

C D|e =(2.56/b)*[r |e *AR*cos\]/[ /+ 

Cfc/Cfj : Account for compressibility effects above Mach * 1 
(Datcom, Figure 4.1.5.1-15). 


Cq j e : Leading edge bluntness term. 

The C Do for the wing is listed in Table ADI. 

BODY: 

SUBSONIC: 

C Do =C Df +C Db 

C Df = 1 .02*C f *( 1 + 1 .5/( 4v*/) 2+ 7/( V*0 4 ) 5 ref /S B 

C Db =0.029*(d b /d) 3 /(C Df ) 0 - 5 

where C Db : Base pressure drag coefficient. 

TRANSONIC: 

^Do _ ^Df + ^Db + ^Dp + *“Dp 
C Dr l.02*C f *[C fc /C f ]* S ref /S B 

C D p- 1 ,02*C f [ 1 ,5 /( Vv> 2 *7/( fc/4 ) 4 J* S ref /S B 
C 0B=- c pb* (d D /d)2 

C pb : Base pressure coefficient (Datcom, Figure 4.2.3.1-17). 

C D ^: Subsonic Pressure Drag. 

C Dp : Supersonic Wave Drag (Datcom, Figure 4.2.3.1-18). 
SUPERSONIC: 

c Do =c f* s ref /5 B +c DN2 +c DA +c DA(NC) + c Db 
C DN2 : Nose wave drag (Datcom, Figure 4.2.3.1-50). 

C DA : Body after body wave drag (Datcom, Figure 4.2.3. 1 -[36-38]). 



CdA(NC)' Interference drag coefficient acting on after body due to 
center body and nose (Datcom, Figure 4.2.3.1-54). 

: Base pressure drag coefficient (Datcom, Figure 4.2.3.1-44) 

The C Do for the body is listed in Table AD2. 

The final C Do value can be obtained from the following equation: 

^Do = ^Do Wing + ^Dotody *^ref 
The (C Do )y ota j is listed in Table AD3 

WING-BODY COMBINATION: 

5-UB5- QN 1& 

C D <Do^*C L 2*K^C L -C Lm |n] 2 

K‘'=1/[iT*AR*e] 

e=e^l 1 -(d/b) 2 ] 

K’l Inviscid drag due to lift (induced drag). 

Viscous drag due to lift due to flow seperation and increased skin 
friction (Nicolai, Figure 1 1.6). 
e : Wing efficiency factor. 

/■. Wing planform efficiency factor (Nicolai, Figure 1 1.5). 

TRANSONIC: 

The values for C L and K have been estimated using Figure AD7 and 
Figure AD8. In these graphs the values of C L and K were estimated 



by simply connecting the subsonic and supersonic values and then just 
reading the corresponding C L and K for the transonic region. 


SUPERSONIC: 


c Do + k * Kl c Lmin^ 

K-1/Cl - AN for subsonic leading edge. 


K= 1 /C 


Loc 


AN=(AN/AN^)*AN^ 0 

AN =1/0.. -(k’+K") 

K: Drag due to lift factor. 

AN: Leading edge suction parameter. 

C|_ : supersonic lift curve slope (Nicolai, Figure 1 1.2). 

(AN/fcN ): Nicolai, Figure 11.7. 

>W‘0 

Tables AD4,AD5, and AD6 list values for C|_ /degree, A N, and 
K, respectively, for the different Mach numbers. 


US 





CALCULATIONS FOR THE SUBSONIC REGION} 
program AAE515H; 
const 
s = 722.5; 
pi = 3.141593; 
cdo = 0.0129385; 
xk = 0.3103235; 


var 

xm, cl, cd, d, v, h, t, r, xx, ala, da, W : real; 
outfile : text; 
tt : boolean; 
aaa : integer; 
oegin 
show text; 

rewrite(outfi1e, 'hani:out4'); 
tt := false; 

while not tt do 
begin 

write! n(*input W,h,t,r'); 

readln(W); 

readln(h); 

readln(t); 

readln(r); 

writelnCw = w : 1 : 2, * lbs’); 
write!n(outfile, *w = ', w : 1 : 2, ’ lbs'); 


xm := 0.1; 

da := 1.6566 * pi / 180.0; 
while (xm <= 0.9) do 


begin 

v := xm * sqrtC 1.4* 1716 * t); 
d := 2 * (w) / r / s / (v * v); 
cd := cdo * 1.1 + (xk * d * cl); 
d := cd * s * (v * v) * r * 0.5; 
xx := d / cd; 


ala := cl / cla; 

Y'/rite(outfile, h : 1 : 1, * xm : 1 : 1, ' v : 1 : 2, ' ’, cl : 1 : 5, 1 ’); 

writeln(outfile, cd : 1 : 5, 1 xx : 1 : 5, ‘ d : 1 : 2, ' ala : 1 : 3); 

write(h : 1 : 1, ‘ xm : 1 : 1, ’ ’, V : 1 : 2, ’ ’, cl : 1 : 5, ' ’); 
writeln(cd : 1 : 5, ’ xx : 1 : 5, ' d : 1 : 2, ' ala : 1 : 3); 
xm := xm + 0.1 ; 


)n 


xm := xm + 0.1; 

end; 

writelnO hit 1 to exit'); 
readln(aaa); 
if (aaa = 1) then 
tt := true; 
end; 

?nd. 


\U 


CALCULATIONS FOR THE TRANSONIC REGIONS} 


irogram AAE515H; 
const 
s = 722.5; 


var 

xm, cl, d, cd, cdo, xk, v, h, t, r, xx, ala, cla, w : real; 
infile, outfile : text; 
tt : boolean; 
aaa : integer; 

procedure detcdoxk (xm : real; 

var cdo, xk, cla : real); 

begin 

if (xm = 9.0) then 
begin 
xk := 0.31; 
cdo := 0.0147191; 
cla := 1.683; 
end; 

if (xm = 10.0) then 
begin 

xk := 0,30; 
cdo := 0,034345; 
cla := 1.734; 

end; 

if (xm = 1 1.0) then 
begin 
Xk := 0.285; 
cdo := 0.039222; 
cla := 1.753; 
end; 

if (xm = 12.0) then 
begin 
xk := 0.28; 

Cdo ;= 0.039493; 
cla := 1.722; 
end; 
end, 

segin 

reset(infile, 'Papers:datar);} 


rewiite(outf11e, 'HaniicutS'); 

It := false; 

while not tt do 
begin 

writelnCinput W,h,t,r'); 

readln(w); 

readln(h); 

readln(t); 

readln(r); 

writeln(outfile, 'w = w : 1:2/ lbs’); 
writelnCw = *, w : 1 : 2/ lbs'); 
xm := 9.0; 

while (xm <= 12.0) do 

begin 

v := xm * sqrtC 1 .4 * t * 1716) / (1 e0 1 ); 
cl := 2 * (w) / r / s / (v * v); 
detcdoxk(xm # cdo, xk, da); 
cd := cdo * 1.1 + (xk * cl * cl); 
d := cd * s * (v * v) * r * 0.5; 
xm := xm / led; 
xx := cl / cd; 

ala := cl / (cla * pi) * 180.0; 

write(outfi1e, h: 1 : 1/ xm:1:1/ v : 1 : 3/ cl : 1 : 5/ '); 

writelnfoutfile, cd : 1 : 5/ *, xx : 1 : 5/ d ; 1 : 2/ ala : 1 : 2); 
{v/ri te(outfile);} 

write(h: 1:1/ xm : 1 : 1/ > : 1 : 3/ *, cl : 1 : 5/ '); 
write1n(cd : 1 : 5/ xx : 1 : 5/ *, d : 1 : 2/ ala : 1 : 2); 

{writeln;} 

xm := xm * 1e01; 
xm := xm + 1.0; 
end; 

writelnOHit 1 to exit'); 
reedln(aaa); 

if aaa = 1 then 
tt := true; 
end; 

2nd. 



ex 


CALCULATIONS FOR THE SUPERSONIC AND HYPERSONIC REIGQNS} 

>rogram AAE515H; 
const 
s = 722.5; 
pi = 3.141593; 


var 

xm, cl, d, cd, cdo, xk, v, xx, ala, da, h, t, r, W : real; 
infile, outfile : text; 
tt : boolean; 
aaa : integer; 

rocedure detcdoxk (xm : real;} 
var cdo, xk : real);} 
begin} 

vl IW,/ 

)egin 

reset(infile, ’Papers:datar);} 
showtext; 

rewri te(outfile, ‘Hani:out6‘); 
tt := false; 

while not tt do 
begin 

writelnCinput W,h,t,r'); 

readln(w); 

readln(h); 

readln(t); 

readln(r); 

write1n(outfi1e, 'w = ', w : 1 : 2, 1 lbs’); 
wntelnCw = w ; 1 : 2, ‘ lbs'); 
xm := 1 5.0; 

while (xm <= 100.0) do 

begin 

v := xm * sqrt( 1 .4 * t * 1 7 1 6) / ( 1 e0 1 ); 
d := 2 * (w) / r / S / (v * v); 
detcdoxk(xm, cdo, xk);} 
if (xm = 15.0) then 
begin 
xk := 0.35; 

Cdo := 1.3 * 0.037375; 


\M 


Cla := 2.2355; 

end; 

if (xm = 20.0) then 
begin 
Xk := 0.430; 
cdo := 1.3 * 0.02S804; 
cla := 2.092; 
end; 

if (xm = 25,0) then 
begin 

xk := 0.5210; 

cdo := 1.3 * 0.025900; 

cla := 1.9193; 

end; 

if (xm = 30.0) then 
begin 
xk := 0.602; 
cdo := 1.3 * 0.02420; 
cla := 1.662; 
end; 

if (xm = 35.0) then 
begin 

Xk := 0.7294; 

Cdo := 1.3 * 0.022705; 
cla := 1.371; 

end; 

if (xm = 40.0) then 

begin 

xk ;= 0.SS02; 

Cdo := 1.3 * 0.21767; 
cla := 1.1361; 

end; 

if (xm = 45.0) then 
begin 
xk := 1.02; 

cdo := 1.3 * 0.02077; 
da := 0,980; 

end; 

if (xm = 50.0) then 
begin 

xk := 1.161; 

cdo := 1.3 * 0.020104; 


end; 

if (xm = 55.0) then 
begin 
xk := 1.294; 
cdo := 1.3 * 0.019427; 
cla := 0.773; 
end; 

if (xm = 60.0) then 
begin 
xk := 1.422; 
cdo := 1.3 * 0.018758; 
cla := 0.703; 
end; 

if (xm = 65.0) then 
begin 

xk := 1.55; 

cdo := 1.3 * 0.0183272; 
cla := 0.645; 

end; 

if (xm = 70.0) then 
begin 
Xk := 1.672; 

cdo := 1.3 * 0.0176583; 
cla := 0.593; 

end; 

if (xm = 75.0) then 
begin 
xk := 1.799; 
cdo := 1.3 * 0.017373; 
Cla := 0.556; 
end; 

if (xm = 80.0) then 
begin 
xk := 1.922; 
cdo := 1.3* 0.017012; 
Cla := 0.5203; 
end; 

if (xm = 85.0) then 
begin 
xk := 2.049; 

cdo := 1.3* 0.0166643; 
cla := 0.488; 

end; 



if (xm = 90.0) then 
begin 
xk := 2.1 S3; 
cdo := 1.3 * 0.01639; 
els := 0.458; 
end; 

if (xm = 95.0) then 


begin 

xk := 2.32; 

cdo := 1.3 * 0.016125; 
cl a := 0.431; 


end; 

if (xm = 100.0) then 
begin 
xk := 2.457; 

Cdo := 1.3 * 0.015517; 

Cla := 0.407; 
end; 

cd := cdo / 1.3 * 1.25 + (xk * cl * cl); 
d := cd * s * (v * v) * r * 0.5; 


xm := xm / (led); 
xx := cl / cd; 

ala := cl / (cla * pi) * 180.0; 

write(outfi1e, h : 1 : 1/ xm:1:1/ *, 7 : 1 : 3 / cl : 1 : 5/ '); 

write!n(outftIe, cd : 1 : 5/ xx : 1 : 5/ d : 1 : 2/ *, ala : 1 : 2); 
write(outfile); 


•vrite(h : 1 : 1/ 
■vn tainted : 1 ; 5, 
n nteln;} 


xm: 1 : 1/ *, v : 1 : 3/ \ cl : 1 : 5/ *);} 
xx : 1:5/ *, d : 1 : 2/ ala : 1 : 2);} 


xm := xm * led; 
xm := xm + 5.0; 


end; 

writelnCHJt ! to exit’); 
readln(aaa); 

if aaa = 1 then 


U := true; 

end; 

nd. 
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program INLETAREA 

d i men s i on H ( 30 ) , D ( 30 ) , T ( 30 ) , MACH ( 30 ) , VEL < 30 ) , BLD ( 30 ) 
d i men s i on HM (30) , DM < 30 ) , TM ( 30 ) , W 1 R < 30 ) , AREA < 30 ) , BARE A (30) 
read MACH,Q 

open (un i t = l , f i 1 e= ' I NLETAREA. DAT ' ) 


ALT =35000. 
do 1=1,13 

ALT =ALT +5000 . 
H ( I ) =ALT 
T (I) =389.99 

end do 

T< 10) =394. 69 
T (1 1 ) =402.48 
T< 12) =410. 64 
T (13) =418. 79 

D < 1 ) =. 00058727 
D (2) =.00046227 
D (3) =. 00036391 
D (4) =.00028652 
D (5) =.00022561 
D (6) =. 00017767 
D (7) =.00013993 
D (8) =.0001 1022 
D (9) =. 00008683 
D< 10) =.00006771 
D < 1 1 ) =. 00005253 
D< 12) =.00004097 
D ( 1 3 ) = . 000032 1 1 


3 — o . 
write(6 
■f or ma t < 

wri te (6 
wri te (6 
wri te (6 
wr i te ( 1 
wr i te d 
wr i te d 
writed 
format ( 

f ormat < 
do 1=1 , 

wr i te (6 


f 1 ow 


, 6X , ' Mass 

y - j 

,7X , ' (f t**2) 


» ' '> 

, 20 ) 

,30) 

,5) 

,5) 

,20) 

, 30) 

,5) 

' Mach # ' , 5X , ' A1 1 i turJe 
, 'Corrected Area',8X, 

14X , ' ( f t ) ' , 8X , ' (1 b/sec ) 

13 

J=J + 1 

, *) 'Enter Mach #, W 1 R , % bleed at 
read (5,*) MACH (J ), W1R ( J ), BLD <J ) 

HM ( J ) =H ( I ) 

DM ( J ) =D ( I ) 

TM ( J ) =T ( I ) 

VEL ( J ) =MACH ( J ) *SQRT ( 1 . 4* 1 7 1 6*TM ( J ) ) 

AREA ( J ) =W 1 R ( J ) / ( 32 2*DM ( J ) * VEL ( J ) ) 

BARE A ( J ) = AREA ( J ) *BLD ( J ) +AREA ( J ) 

Y=BAREA ( J ) n . 

40 ) MACH ( J ) , HM ( J ) , W 1 R ( J ) , AREA ( J ) , BAREA ( J ) 


hi- ea 




,7 X , ' (ft-* *2) 


,9X , 


(ft#*! 


, H ( I ) 


ft 


wri te <6 

writed, 40 ) MACH ( J ) , HM ( J ) , W 1 R < J ) , AREA ( J ) , BAREA ( J ) , Y 
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40 


-format (1 X , F4 ■ 2 , SX , FS. 1 , 5X , F6. 1 , 1 OX , F6. 3 , 8X , F6. 3 , 10X , F5. 3) 
write <6, *) 'Enter 1 to enter another Mach# at same altitude else 
read (5,*) Q 

if <Q. eq. 1 ) goto 10 

end do 
cl ose < 1 ) 
stop 
end 


( 





program INLETANGLEM6 

dimension M(15) ,MN(15> , BET (15) , THET (15) , XT (15) , FOR (15) 
real M,MN 

open (uni t=l , f i 1 e=’ INLANGM6. DAT’ ) 
write (6,*) ’Enter Theta, Ml, NR’ 
read (5,*) THETA, M ( 1 ), NR 
M ( 1 ) =6. 

NR=8 

THET ( 1 ) =THETA*3. 141592654/180. 

BET ( 1 ) =BETA*3. 141592654/180. 

TT=TAN (THET ( 1 ) ) /2. 


do 30 1=1, NR 

write(6,*) ’Enter Theta -for Ml=6 and NR=8’ , I 
read (5,*) THETA 
BET A=THET A+3 . 

THET ( I ) =THETA*3. 141592654/180. 

BET ( I ) =BETA*3. 141592654/180. 

TT =T AN ( THET ( I ) ) /2. 

TB=0. 

i-f (ABS(TT-TB) . It. .00001) goto 20 
BET ( I ) =BET ( I ) + (TT-TB) *1.7 

A=(M(I)**2.*(1.4+C0S<2.*BET<I) ) )+2. >*TAN(BET(I> ) 
TB=(M(I) **2. *SIN(BET(I) ) **2. -1 . ) /A 
goto 10 

XT ( I ) =4. /TAN (BET (I) ) 

MN(I)=M(I) *SIN(BET(I) ) 

MN ( 1+1 ) =SQRT ( (MN ( I ) ##2. +2. / . 4) / (2. 8/. 4*MN ( I ) **2. -1 ) ) 
M ( 1+1 ) =MN ( 1+1 ) / ( S I N ( BET ( I ) -THET ( I ) ) ) 

B= ( 1 . 1 6666666667 *MN ( I ) **2. -. 166666666667) **-2.5 
POR ( I ) = ( 1 . 2*MN ( I ) **2. / (1 . +. 2*MN ( I ) **2. ) ) **3.5*B 
THET ( I + 1 ) =THET ( I )- 
BET ( I + 1 ) =THET ( I + 1 ) 
continue 


wr i te (6, 50) 

■format (’ THETA’ ,4X, ’BETA’ , 7X , ’ M ( I ) ’ , 7X , ’M(I + 1) ’ ,5X, ’FOR’ , 8X, ’ XT’ ) 
wri te ( 1 , 50) 

TPOR-1. 
do 40 1=1, NR 

BETA=BET('l) /3. 141592654*180. 

THETA=THET( I) /3. 141592654*180. 

wri te (6, 60) THETA, BETA, M ( I ) ,M(I+1) , POR < I ) ,XT(I) 
format ( 1 X , F5. 2, 5X , F5. 2, 5X , F6. 3, 5X , F6. 3, 5X , F6. 5, 5X , F5. 2) 
wri ted, 60) THETA, BETA, M ( I ) ,M(I + 1) , POR ( I ) , XT (I) 

TPOR-POR ( I ) *TP0R 
continue 

C= ( 1 . 16666666667*M (NR+1 ) **2. -. 166666666667) **-2.5 
TPTR=( 1 . 2*M (NR+1 ) **2. / ( 1 . +. 2*M (NR+1 ) **2. ) ) **3. 5*C*TP0R 
wr i te (6, 70) 

format!/,’ TPOR’ , 7X , ’ TPTR’ ) 
write (6, 80) TPOR, TPTR 
format ( 1 X, F6. 5, 6X , F6. 5) 

wri te (6, *) ’ 

wr i te (1 , 70) 

wri te(l, 80) TPOR, TPTR 

wri ted, *) ’ 

goto 1 
close ( 1 ) 
stop 

end \'L% 








THETA 

BETA 

MCI) 

M < I + 1 > 

POR 

XT 

5 - 00 

13. 16 

6 . 000 

5.315 

„ 96624 

17.11 

6 „ 00 

15.20 

5. 315 

4.647 

. 95987 

14.73 

7 . 25 

17.79 

4.647 

3. 997 

. 95305 

12.46 

a. 75 

21.13 

3.997 

c* m y 2 

. 94735 

10. 35 

1 0 . 25 

25.28 

3.372 

2.790 

. 94763 

8.47 

1 2 . 00 

30.92 

2. 790 

2. 240 

. 94947 

6 . 68 

13.25 

38 . 60 

2 . 240 

1 . 730 

. 95887 

5.01 

. n C - 

J> X. ■ ^1. V_J 

49.40 

1 . 730 

1 . 290 

. 97690 

3. 43 

TP OR 

TPTR 





. 70576 

.69241 





THETA 

BETA 

M < I ) 

M < I + 1 ) 

POR 

XT 

5.25 

13.36 

6.000 

5.282 

.96146 

•16.84 

6.25 

15.47 

5.282 

4.594 

. 95595 

14. 45 

7.25 

17.94 

4.594 

3. 956 

. 95443 

12. 36 

a . 50 

21.07 

3.956 

3. 35 7 

. 95252 

10.38 

9.75 

24.91 

T ?C! -7 
*-.* m / 

2. 806 

. 95462 

8.61 

1 1 . 00 

29.82 

2.806 

2. 299 

. 95932 

6.98 

12.25 

36 . 63 

2.299 

1.822 

. 96497 

5. 38 

13.75 

48.54 

1 . 822 

1 . 323 

' . 96644 

*-> « 5 o 

TF'OR 

TPTR 





. 71362 

.69591 





THETA 

BETA 

M ( I ) 

M < I + 1 > 

POR " 

______ 

5 - 25 

il O* m O G> 

6.000 

5. 282 

.96146 

16.84 

6.25 

15. 47 

5.282 

4.594 

. 95595 

14.45 

7.25 

17.94 

4 . 594 

3. 956 

. 95443 

12. 36 

8 . 50 

2 1 . 07 

3.956 

3. 357 

.95252 

10. 38 

9.75 

24.91 

3.357 

2. 806 

. 95462 

8.61 

1 1 . 00 

29.82 

2.806 

2.299 

. 95932 

6.98 

13.75 

38.25 

2.299 

1 . 761 

. 95222 

5.07 

12.25 

48.32 

1 . 761 

1 . 322 

. 97668 

3 . 56 

TP OR 

TPTR 





.711 65 

. 694 1 0 





THETA 

BETA 

M ( I ) 

M < I + 1 > 

POR 

Xj 

5 - 25 

1 3. 36 

6. 000 

5.282 

.96146 

1 6 .84 

6 . 25 

15.47 

5.282 

4.594 

. 95595 

14.45 

7 - 25 

1 7 . 94 

4.594 

3. 956 

. 95443 

12. 36 

S. 50 

21 . 07 

3.956 

3.357 

. 95252 

1 0 . 38 

9.75 

24. 91 

•■1"* m v>iJ / 

2. 806 

.95462 

8.61 

1 1 . 00 

29.82 

2 . 806 

2- 299 

. 95932 

6.98 

1 3 . 00 

37.43 

2 • P99 

1.791 

. 95888 


1 3 . 00 

48.41 

1.791 

1 “TOT 

J. m *.«• Xm 

.97189 

3 , 55 

TP OR 

TPTR 





.71311 

. 69540 





\ %<H 



THETA 

BETA 

M (I) 

M < .1 + 1 ) 

FOR 

XT 

5.00 

13. .16 

6.000 

5.315 

. 96624 

.17.1 1 

5 . 00 

14. 39 

5.315 

4.756 

. 97563 

15.59 

5 . 00 

15. 67 

4.756 

4.288 

. 98195 

14.26 

5 . 00 

1 7 . 02 

4.288 

3. 888 

. 98635 

13. 07 

5 . 00 

18. 45 

3.888 

3. 540 

. 98950 

11.99 

8.33 

22. 72 

3.540 

3.035 

. 96600 

9.55 

8 . 33 

25. 65 

3.035 

2.617 

.97687 

8.33 

S m v>s!> 

29. 10 

2.617 

m X~\-t / 

. 98384 

7. 19 

TPOR 

TPTR 





. 83877 

. 505 1 6 





THETA 

BETA 

M< I ) 

Hd + 1) 

FOR 

. - t - 

6 . 00 

13.98 

6 . 000 

5. 182 

.94491 

16.06 

6 . 00 

15.47 

5. 182 

4.541 

. 96247 

14.45 

6 . 00 

17. 05 

4.541 

4 . 020 

. 97355 

1 3 . 04 

6.00 

18. 74 

4.020 

3. 585 

.98085 

1 1 . 79 

6 . 00 

20.56 

3.585 

3.213 

. 98583 

1 0 . 66 

1 0 . 00 

25. 97 

3.213 

2.677 

. 95641 

8.21 

1 0 . 00 

30. 03 

2.677 

D 077 
X. • X- N-‘ / 

. 97180 

6. 92 

:t 0 . 00 

35 * 

— r-y-ry 

■»_ M Jb. •_* / 

1.856 

.98102 

5 « 67 

TPOR 

TPTR 





. 78061 

. 61478 





THETA 

BETA 

M ( I > 

M < I + 1 ) 

POR 

XT 

7 . 00 

14. 84 

6 . 000 

5.048 

.91779 

15. 10 

7 „ 00 

16. 60 

5.048 

4.334 

. 94705 

13.42 

7 . 00 

18. 49 

4.334 

3.769 

. 96435 

1 1 . 96 

7 . 00 

20.55 

3.769 

3. 306 

. 97519 

1 0 . 67 

7 . 00 

22.82 

3.306 

2. 916 

. 98231 

9.51 

11. 67 

29.56 

2.916 

2.356 

. 9480.1 

7.05 

11.67 

35.24 

2.356 

1.894 

.96771 

5. 66 

11.67 

43. 62 

1 . 894 

1 . 478 

.97815 

4 . .20 

TPOR 

TPTR 





. 72053 

. 67477 





THETA 

BETA 

M < I ) 

M <1+1 ) 

POR 

XT 

3 . 00 

15.73 

6 „ 000 

4. 915 

. 885.28 

14.21 

8 „ 00 

17.77 

4.915 

4. 133 

. 92986 

12.48 

8 . 00 

19.99 

4. 133 

‘TT t cr t; T s 

. 95484 

10.99 

8. 00 

22.46 

* i— J •— 1 

3. 050 

. 96972 

9 . 67 

8 . 00 

25.27 

3 . 050 

2. 645 

. 97909 

8 . 48 

1 •-» u 3 o 

33.59 

2.645, 

2. 063 

. 94092 

6.02 

1 s!> » •«' 

4 1 . 84 

2 . 063 

1.571 

„ 96390 

4.47 

i 3 n 33 

60. 88 

1.571 

1.019 

. 96474 



TPOR TPTR 

. 65296 . 65296 



THETA 

BETA 

M<I) 

MU + 1) 

POR 

XT 

7.00 

14. 84 

6 . 000 

5. 048 

.91779 

15. 10 

7 , 00 

16. 60 

5.048 

4.334 

. 94705 

13.42 

7 „ 00 

18. 49 

4.334 

3. 769 

. 96435 

1 1 . 96 

7 . 00 

20. 55 

3. 769 

3 . 306 

. 97519 

10.67 

7 00 

22. 82 

3 u 306 

2.916 

. 9823 1 

9.5.1 

10. 00 

28.01 

2. 916 

2 * 4c:'d 

„ 96550 

7 . 52 

1 1 . 00 

33 . 56 

2.436 

1 . 990 

. 97045 

6 . 03 

1 4 . 00 

44. 26 

1 . 990 

1.477 

. 96094 

4.10 

TP OR 

TPTR 





. 72295 

. 67730 





THETA 

BETA 

MU) 

M < 1 + 1 ) 

POR ’ 

XT 

5 . 00 

13. 16 

6.000 

5. 315 

. 96624 

17.11 

6 . 00 

15.20 

5.315 

4. 647 

. 95987 

1 4 . 73 

7 . 00 

17.58 

4 . 647 

4. 019 

. 95725 

12.62 

3 . 00 

20. 40 

4.019 

3. 443 

. 95793 

10. 76 

9 . 00 

23.79 

3.443 

2.917 

. 96100 

9. 07 

1 0 . 00 

28. 00 

•2.917 

2.437 

. 96547 

7.52 

1 1 . 50 

34.04 

2. 437 

1 . 970 

. 96663 

5. 92 

1 3 . 50 

44.09 

1 . 970 

1 . 479 

. 96520 

4. 13 

TF'OR 

TPTR 





. 73621 

. 68940 





THETA 

BET A 

Mil) 

M U + l ) 

POR 

XT 

5 . 00 

13. 16 

6. 000 

5. 315 

. 96624 

17.11 

5 . 75 

14.99 

5.315 

4. 674 

. 96424 

14.94 

6. 50 

17.09 

4.674 

4.085 

. 96446 

13.01 

7 . 50 

19.74 

4. 085 

3 . 531 

.96313 

11.15 

8.75 

23. 12 

3. 531 

3. 003 

.96147 

9.37 

1 0 . 00 

27.36 

3.003 

2.507 

.96299 

7 . 73 

10.75 

32.49 

2.507 

2. 061 

. 97042 

6.28 

11.25 

39. 49 

2. 061 

1 . 650 

. 97748 

4 . 85 

TPOR 

TPTR 





. 76009 

. 66590 





THETA 

BETA 

MU) 

MU + 1) 

POR 

XT 

4 . 00 

12.37 

6 . 000 

5.450 

.98180 

1 8 . 24 

6 . 00 

14. 93 

5. 450 

4.754 

. 95715 

1 5. 00 

3 . 00 

18. 16 

4. 754 

4.012 

. 93553 

12.1 9 

1 0 . 00 

22. 19 

4. 012 

3.295 

. 92489 

9.81 

1 2 . 00 

27. 34 

3.295 

2 . 633 

. 92552 

7 . 74 

1 0 . 00 

3.0. 46 

2 . 633 

2. 199 

. 97287 

6 . SO 

1 3 . 00 

38. 97 

2.199 

1 . 706 

. 96226 

4 . 94 

1 7 . 00 

62.94 

1 . 706 

. 966 

. 92369 

2.04 

TPOR 

TPTR 





. 65075 

. 65078 





THETA 

BETA 

M < I ) 

M ( l' + l ) 

POR 

XT 

5 00 

13.16 

6 . 000 

5. 315 

„ 96624 

17.11 

6 . 00 

15.20 

5. 315 

4.647 

. 95987 

1 4 „ 73 

7 „ 25 

1 7 . 79 

4 . 647 

3. 997 

. 95305 

1 2 . 46 

8 . 75 

2 1 1 3 

3. 997 

3. 372 

94735 

1 0 . 35 

1 0 . 25 

25.28 

■i *r ^ 5* 

2. 790 

. 94763 

8. 47 

1 3 . 00 

3 1 . 90 

2.790 

2. 193 

. 93773 

6 . 43 

1 2 00 

37 . 98 

2. 193 

1 „ 739 

. 96983 

5. 12 

12.75 

49.86 

1 „ 739 

1 . 278 

. 97389 

3 a 3 / 


TP OR TPTR 

. 70.283 „ 69033 


\1 I 



THETA 

BETA 

MCI) 

M <1 + 1 ) 

POR 

XT 

4 . 00 

12.37 

6. GOG 

5 . 450 

. 98180 

18.24 

5 . GO 

14.12 

5 . 450 

4. 867 

. 97393 

15.90 

6 . GO 

16. 19 

4.867 

4.287 

. 968 1 9 

13. 78 

7 - GO 

18.64 

4.287 

3. 731 

. 96534 

11. 86 

8 . GO 

21.55 

o m 7 •«' 1 

3.211 

. 96521 

10. 13 

8 . GO 


3.21 1 

2.782 

. 97622 

8 . 89 

1 0 . GO 

29.10 

2.782 

2.325 

.96917 

7. 19 

12. GO 

36. 01 

2.325 

1 . 854 

. 96611 

5.50 

TPOR 

TPTR 





. 78848 

.62173 





THETA 

BETA 

M < I ) 

M < 1 + 1 ) 

POR 

XT 

5. GO 

1 0 • i to 

6 . 000 

5. 315 

. 96624 

17.11 

6 . GO 

15.20 

5.315 

4. 647 

. 95987 

14.73 

7 . GO 

17. 58 

4.647 

4.019 

. 95725 

12. 62 

8 . GO 

20 . 40 

4.019 

o ■ *4 

. 95793 

10.76 

9 . GO 

23. 79 

3. 443 

2.917 

. 96 1 00 

9 ,, 07 

9 GO 

27. 10 

2.917 

2.485 

. 97407 

7.82 

12. GO 

33.98 

2.485 

1 . 989 

. 96098 

5. 94 

14. GO 

44.27 

1.989 

1 . 477 

. 96095 

4 . 1 0 

TPOR 

TPTR 





. 73517 

. 68881 





THETA 

BETA 

MCI) 

M < I + 1 ) 

POR 


5 . GO 

13. 16 

6 . 000 

5.315 

. 96624 

1 7.11 

5 . GO 

14.39 

5.315 

4.756 

. 97563 

15.59 

6 . GO 

16. 47 

4.756 

4. 196 

. 97008 

1 3 . 53 

6. GO 

18. 12 

4. 196 

3. 733 

. 97854 

12.23 

7 . GO 

20.70 

3 » 733 

3.276 

. 97580 

10.58 

7. GO 

22.99 

3.276 

2.891 

. 9827 1 

9. 43 

1 0 . GO 

23.20 

2.891 

2. 415 

. 96622 

7 . 46 

12.00 

34.82 

2.415 

1.931 

. 96329 

b u / 5 

TPOR 

TPTR 





. 79869 

. 60 1 50 





THETA 

BETA 

MU) 

M <1 + 1 ) 

""for"’ ' 

XT 

4 . GO 

12. 37 

6.000 

5. 450 

. 98180 

18.24 

5„ 00 

14. 12 

5 . 450 

4. 867 

.. 97393 

15.90 

7 . GO 

1 7 . 02 

4 . 867 

4. 193 

. 95183 

13.07 

8 . GO 

19.79 

4. 193 

3. 580 

. 95315 

1 1 . 1 2 

i G . GO 

23. 99 

3 . 580 

2.967 

. 94318 

8. 99 

1 0 . GO 

27.62 

2.967 

2. 478 

. 96405 

7.64 

1 1 . GO 

33. 07 

2.478 ' 

2.025 

. 96933 

6.14 

1 3 . GO 

42.23 

2.025 

1 . 550 

. 96733 

4.41 

TPOR 

TPTR 





. 73964 

. 67549 











THETA 

BETA 

Mil) 

M i I + 1 ) 

FOR 

XT 

4.50 

12.76 

6.000 

5. 383 

. 97473 

17.66 

5 . 50 

1 4 . 65 

5. 383 

4. 756 

.96723 

1 5 . 30 

6 . 50 

1 6 . 88 

4 . 756 

4. 151 

. 96281 

13. 18 

7 . 50 

19.51 

4. 151 

3. 584 

.96159 

11.29 

8 . 50 

22. 65 

3.5B4 

3. 061 

. 96301 

9.59 

9.50 

26.51 

a 0 U) .1 

2. 580 

. 96625 

8 . 02 

1 0 . 50 

3 1 . 47 

2 . 580 

2. 133 

. 97044 

6.53 

12. 50 

39. 56 

2. 133 

1 . 666 

. 96797 

4.84 

TPOR 

TPTR 





. 76295 

. 66337 

■ 




THETA 

BETA 

Mil) 

Mi 1 + 1 ) 

POR 

XT 

4 . 50 

.12. 76 

6 . 000 

5.383 

. 97473 

1 7 . 66 

5 . 50 

14. 65 

5 . 383 

4. 756 

. 96723 

1 5 . 30 

7 . 00 

17. 30 

4.756 

4. 105 

. 95462 

12. 84 

8.50 

20. 52 

4. 105 

3. 474 

.94785 

1 0 . 68 

1 1 . 00 

25. 43 

3. 474 

2.825 

. 93239 

8.41 

1 1 . 50 

30 . 1 3 

2.825 

2.291 

. 95349 

6 . 89 

12.00 

36.48 

2.291 

1.825 

. 96709 

5.41 

1 3 . 00 

47.32 

1 . 825 

1 „ 358 

.97145 

3 .69 

TPOR 

TPTR 





. 71250 

. 68975 





THETA 

BETA 

Mil) 

M(I+1) 

POR 

XT 

4 . 50 

12.76 

6 . 000 

5. 383 

. 97473 

17.66 

5. 00 

14.25 

5. 383 

4.812 

. 97479 

15. 75 

5.75 

16. 13 

4.812 

4.265 

. 97255 

13. 84 

6 . 75 

IS. 50 

4.265 

O a 7 

. 96901 

11.95 

8 . 00 

21.54 

3 . 733 

3.212 

.96517 

10 , 13 

9.25 

25 a vlj O 

3.212 

2.717 

. 96464 

8.46 

9.50 

29.21 

2.717 

2. 293 

. 97464 

7. 15 

11.25 

35. 69 

2 *7 9 3 

1 . 356 

. 97237 

5.57 

TPOR 

TPTR 





. 79010 

. 62228 





THETA 

BETA 

Mil) 

Mi 1 + 1) 

POR 

Xf 

4 . 50 

12.76 

6 „ 000 

5. 383 

. 97473 

1 7 . 66 

6 . 00 

1 5 . 06 

5.383 

4. 700 

. 95853 

1 4 . 8 6 

7 . 50 

17.86 

4.700 

4. 017 

. 947 1 1 

12. 41 

9 . 50 

21 . 72 

4.017 

3. 334 

. 93396 

1 0 . 04 

1 1 . 00 

26. 18 

sl‘« a 334 

2.718 

. 93874 

8.14 

12.00 

31.56 

2.718 

2 . 182 

. 95236 

6.51 

1 2 . 50 

38. 71 

2.182 

1.71 0 

. 96659 

4 . 99 

1 4 . 00 

53 . 30 

1.710 

1 . 188 

. 965 1 5 

2 . 98 

TPOR 

TPTR 





. 68928 

. 68509 





THETA 

BETA 

Mil) 

M i .1 + 1 ) 

POR 

■ XT 

4 . 50 

12 . 76 

6 . 000 

5. 383 

. 97473 

1 7 . 6 6 

5. 00 

14.25 

5.383 

4. 812 

. 97479 

15. 75 

6 . OO 

16,. 33 

4.812 

4. 242 

. 96914 

1 3 . 65 

8 . OO 

1 9 . 63 

4.242 

3 . 6 .1. 8 

. 95175 

11.21 

1 O . OO 

2 S „ 8 1 

3. 618 

2. 996 

. 94169 

9 . 06 

1 O . OO 

27.41 

2. 996 

2. 502 

. 963 1 9 

7.71 

i l . 00 

32. 79 

2 . 502 

2 . 046 

. 96867 

6,21 

12.50 

41.21 

2 . 046 

. 1.588 

. . 9702 1 

4.57 

TPOR 

TPTR 





. 74709 

. 67212 






\2"b 


THETA 

BETA ■ 

M < I ) 

M < I + 1 ) 

FOR 

XT 

5 . 00 

13.16 

6. 000 

5. 315 

. 96624 

.17.11 

6. 00 

15.20 

5.315 

4 . 647 

. 95987 

1 4 . 73 

7.25 

17. 79 

4 . 647 

3. 997 

. 95305 

12. 46 

S. 75 

21 . 13 

3.997 

3. 372 

. 94735 

1 O'. 35 

10.25 

25.28 

3.372 

2. 790 

. 94763 

8.47 

1 1 . 00 

29.96 

2.790 

2.286 

. 95987 

6.94 

.1 2 . 00 

36. 56 

2. 286 

1 . 820 

. 96726 

5. 39 

14.25 

49.37 

1 . 820 

1 . 300 

. 96275 

3.43 

TF'OR 

TPTR 





. 70931 

. 69469 

- 




THETA 

BETA 

M < I ) 

M(I+1) 

POR 

XT 

5. 50 

13.57 

6. 000 

5. 250 

. 95633 

16.58 

6.50 

15.74 

5 . 250 

4.542 

.95138 

14.19 

S . 25 

18.94 

4. 542 

3. 830 

. 93774 

11. 66 

9. 75 

22.67 

3 . 830 

3. 175 

.93721 

9.58 

11.25 

27.36 

3. 175 

2.581 

. 94224 

7 . 73 

12. 75 

33 - 66 

2.581 

2 . 037 

. 95020 

6 . 0 1 

1 4 . 00 

43.20 

2 . 037 

1.521 

. 95954 

4 . 26 

14. 50 

ft- -X* 

1.521 

2.059 

■X* 

--7. 26 

TF'OR 

TPTR 





.71993 

. 49905 





THETA 

BETA 

M (I ) 

M < I + 1 ) 

POR 

~XT 

5. 50 

13.57 

6 . 000 

5.250 

. 95633 

1 6 . 58 

6 . 50 

15.74 

5.250' 

4. 542 

.95188 

14.19 

7.50 

18.29 

4.542 

3. 894 

.95157 

12.1 0 

a „ 50 

21.31 

3.894 

3. 309 

.95437 

10.26 

9 . 50 

24.96 

3 . 309 

2. 781 

.95915 

8 . 59 

1 0 . 50 

29.56 

2. 781 

2.301 

. 96483 

7.05 

12.50 

36 . 86 

2.301 

1.814 

.96292 

o . 

1 4 . 50 

49.99 

1.814 

1 . 282 

. 96087 

3 . 36 

TP OR 

TPTR 





. 70784 

. 69535 





THETA 

BETA 

M (I) 

M < I + 1 > 

POR 

XT 

5 . 00 

13. 16 

6 . 000 

5.315 

. 96624 

17. 1 1 

6.50 

15.61 

5.315 

4 „ 593 

. 95030 

14.32 

7 . 50 

18. 15 

4 . 593 

3. 934 

. 950 1 6 

12. 20 

8 „ 50 

21.15 

3.934 

w* » 

. 95318 

1 0 . 34 

9 . 50 

24.78 

3. 340 

2. 806 

. 958 1 9 

8 . 66 

10.25 

29 . 1 2 

2 « 806 

2. 334 

. 96633 

7. 18 

1 2 . 25 

36. 14 

2.334 

1 . 851 

. 96388 

5. 48 

1 4 . 50 

48.73 

1 . 851 

i “TOO 

. 9603 1 

3.51 

TP OR 

TPTR 





.71273 

. 69523 





THET A 

BETA 

ii •: i > 

M (I + 1 ) 

POR 

XT 

6 . 00 

1 3 . 98 

6. 000 

5 . 1 82 

. 94491 

16. 06 

7 . OO 

1 6 . 3 .1 

5. 182 

4.437 

. 94336 

1 3 . 67 

■8 . 00 

19.03 

4. 437 

3. 770 

.. 94589 

1 1 . 60 

9 . 00 

22.25 

3.770 

3. 177 

.95104 

9 . 78 

1 0 00 
1 2 „ 00 
1 3 . 00 

26.19 

3. 177 

2„ 649 

»* u / O / 

8.13 

40.26 

2. 649 
2. 125 

2. 125 
1 . 640 

. 95507 
. 96456 

6 . 35 
4 . 72 

1 5 . 00 

60 49 

1 . 640 

1 . 022 

. 95115 

2 . 26 

TPOR 

TPTR 





. 67280 

. 67279 






V-i 4- 



program INLETANGLE 

dimension M(15> ,MN(15> , BET (15) ,THET(15) , XT (15) , FOR (15) 
dimension PR ( 15) , TR ( 15) , DR ( 15) 
real M,MN,MNT 

open <uni t=l , f i le=’ INLANG2.DAT’ ) 
write(6,*) ’Enter Ml, NR’ 
read (5, *) M ( 1 ) , NR 

do 30 1=1, NR 

write(6,*) ’Enter Theta -for Ml and NR’ , I 
read (5, * ) THETA 
BET A=THET A+3 . 

THET ( I ) =THETA*3. 141592654/180. 

BET ( I ) =BETA*3. 141592654/180. 

TT=TAN (THET ( I ) ) /2. 

TB=0. 

if (ABS (TT-TB) . It . . 00001 ) goto 20 
BET ( I ) =BET ( I ) + (TT-TB) *1.7 

A= (M ( I ) **2. * (1. 4+C0S(2. *BET ( I ) ) )+2. ) *TAN (BET ( I ) ) 

TB=(M( I ) **2. *SIN(BET (I) ) **2.-1 . ) /A 
goto 10 

XT ( I )=3. 5/TAN (BET (I) ) 

MN ( I ) =M ( I ) *SIN (BET ( I ) ) 

MN ( 1 + 1 ) =SQRT ( (MN(I) **2.+2. /.4) / (2. 8/ . 4*MN ( I ) **2.-1) ) 

M ( I+l ) =MN ( I+l ) / (SIN (BET ( I ) -THET ( I ) ) ) 

B= ( 1 . 1 6666666667*MN ( I ) **2. — . 166666666667) **-2.5 
PR ( I ) = 1 . 16666666667*MN ( I ) **2. — . 166666666667 
DR ( I ) = (2. 4*MN ( I ) * *2 . ) / (. 4* MN ( I ) **2. +2. ) 

TR(I)=PR(I)/DR(I) 

FOR ( I ) = ( 1 . 2*MN ( I ) **2. / ( 1 . +. 2*MN ( I ) **2. ) ) **3. 5*B 
continue 

wr i te (6, 12) M(l) 
wr i te ( 1 , 12) M(l) 
format (30X, ’Mach # = ’,F5.2) 
wr i t e ( 6 , 15) 

format ( ’ ’ 

, •’ ’ ) 

wri te (6, 50) 
wr i te (6, 15) 

format ( ’ THETA’ ,3X, ’ BETA’ , 6X , ’ M ( I ) ’ , 6X , ’ M ( I+l ) ’ , 4X , ’POR’ ,7X, ’PR’ 

,9X, ’TR’ ,9X, ’ XT’ ) 
wr i te ( 1 , 15) 
wri te ( 1 , 50) 
wr i te ( 1 , 15) 

TP0R= 1 . 

TPR=1. 

TTR= 1 . 

TDR=1 . 

do 40 1=1, NR 

BETA=BET ( I ) /3. 141592654*180. 

THETA=THET ( I ) / 3 . 141592654*180. 

write (6, 60) THETA, BETA, M ( I ) , M ( I + 1 ) , POR ( I ) , PR ( I ) , TR ( I ) , XT ( I ) 

format ( 1 X , F5 . 2 , 4X , F5 . 2 , 4X , F6 . 3 , 4X , F6 . 3 , 4X , F6 . 5 , 4X , F7 . 4 , 4X , F7 . 4 , 4X , F5 . 

wr i te ( 1 , 60 ) THETA , BETA , M ( I ) , M ( I + 1 ) , POR ( I ) , PR ( I ) , TR ( I ) , XT ( I ) 

TP0R=P0R ( I ) *TP0R 
TPR=PR ( I ) *TPR 
TDR=DR (I) *TDR 
TTR=TR ( I ) *TTR 

conti nus - 



C= ( 1 . 16666666667*M (NR+1 ) **2. 166666666667) **-2. 5 
MNT=SQRT ( (MCNR+l) **2. +2. /.4> / <2. 8/ . 4*M (NR+1 ) **2. -1 ) ) 

TPOTR= < 1 . 2*M (NR+1 ) **2. / ( 1 . +. 2*M (NR+l ) *#2. ) ) *#3. 5*C#TP0R 
TPTR= < 1 . 16666666667*M(NR+1 ) **2. 166666666667) *TPR 
TDTR=(2.4*M(NR+1> #*2. ) / ( . 4* M (NR+1 ) **2. +2. > *TDR 
TTTR=TPTR/TDTR 
wri te (6, 17) 

format (’ — **^****^^^S:****^^**** 

? ^ ^ ^ ^ ^ ^ ^ ^ Ur Ur Ur U* Ur Ur Ur Ur ? V 

wri te <6, 70) 

-format (’ TPOR ’ , 6X, ’TPR’ , 8X, ’ TTR ? , 11X, ’ TPOTR’ , 7X, ’ TPTR’ , 8X, ’ TTTR’ 
wri te (6, 80) TPOR, TPR, TTR, TPOTR, TPTR, TTTR, MNT 

format < 1 X , F6. 5, 3X , F8. 4, 4X, F7. 4, 7X, F6. 5, 6X, F8. 4, 4X, F7. 4, 4X , F6. 4) 
wri te ( 1 , 17) 
wri te ( 1 , 70) 

wri te ( 1 , 80) TPOR, TPR, TTR, TPOTR, TPTR, TTTR, MNT 

wr i te (6 , # ) 

wri te ( 1 , *) 

goto 1 

close ( 1 ) 

stop 

end 


FOR MACH #'S M=5.0 & M=6.0 


program INLETLENGTH 
dimension M , BET , THET 

real Ml , MN1 , M2 , MN2 , M3 , MN3 , M4 , MN4 , M5 , MN5 , M6 , MN6 , M7 , MN7 , MWN2 
rea 1 MNN3 , MNN4 , MNN5 , MNN6 , MNN7 
XP I =3. 141593/180. 

write<6,*) ' ENTER THE VALUES FOR THETA1 AND Ml' 
read (5,#) THETA, Ml 

write ( & , * ) 'ENTER THE VALUE OF V AND L1-L3 ' 

read ( 5 , * ) Y , XL1 , XL2 , XL3 , XL4 , XL5 

Ml=5 . 0 

THETA=5 . 25 

Y=3 . 5 

XL1=7 . 270 

XL2=3. 08523 

XL3=1. 7036 117 

XL4=1 . 10302 

XL5=1 .6 

THETA=S . 0 


THESE ARE VARYING LENGTHES 

XL 1=7. 081 003 
XL2=2. 04934137 
XL3=1 .09605 
XL4=0 . 714048 
Pl = l .0 

BETA=THETA+3. 

THET1=THETA*XPI 
BET1=BETA*XPI 
TT=TAN ( THET 1 ) /2 . 

TB=0 . 

if (ABS(TT-TB) .It. .00001) goto 20 
BET1=BET1+(TT-TB) #1 .7 

A=(M‘1**2*(1.4+C0S(2.*BET1) ) +2 . ) *T AN ( BET 1 ) 

XXX=(SIN( BET 1 ) )**2 

TB=<M1**2* ( (SIN(BET1 ) ) **2 )-l . ) /A 

goto 10 


write ( 6 , * ) ' BET 1 =',BET1 

MN1=M1*SIN(BET1 ) 

write ( 6 , * ) ' MN1 = ',MN1 

MN2=SQRT( (MNl*#2+5)/(7*MNl##2-l) ) 

write ( 6 , * ) ' MN2 =',MN2 

write ( 6 , * ) BET 1 , THET 1 

M2=MN2/ ( SI N ( BET i— THET 1 ) ) 

P2=Pl*<l+(2.8/2.4)*(MNl**2-l) ) 

P01=F'l* ( ( 1+0.2*M1**2) ) **3. 5 

P02=P2*( <1+0.2*M2**2) )*'*3.5 

write (6,*) M2 

continue 


ORIGINAL PAGE IS 
OF POOR QUALITY 


do 40 1=1, NR 

BETA=BET( I ) /3. 141592654*180. 

THET A=THET( I )/3. 141592654*180. 
write (6,* ) BETA , THETA , XT ( I ) , M < I ) , M ( I +1 ) 
continue 
BETA^BETl/XPI 

write ( 6 , * ) ' BETA = ' , BET A 


m 



X=Y/TAN ( BETA ) 

Xl = X!_i*COS(THETi) 

Y1=XL1*SIN(THET1 ) 

BET A2=AT AN ( ( Y-Yl ) / ( X-X 1 ) )— THET1 
BETA2=BETA2/XP I 
write (6, 50) X,Y 

format ( iox , * X = * , F6 . 3 , 3X , * Y = *,F6.3) 
write (6 ,60) BETA2 , XI , Y1 

format (IOX , * BETA2= ' ,F6.2,3X , ' Xl=' ,F6.3,3X, ' Yl=' , F6 . 3 ) 

BETA2=BETA2*XF'I 

read ( 5 , * ) 

Al=( ( 2/TAN (BETA2)*(M2**2.*(SIN(BETA2) )**2-l. ) ) ) 
THET2=ATAN(A1/(M2**2.*(1.4+C0S(2*BETA2) ) +2 ) ) 

MNN2=M2*SI N ( BETA2 ) 

write (6,*)' MNN2 = ' , MNN2 

MN3=SQRT ( (MNN2**2+5)/(7*MNN2**2-l) ) 

M3=MN3/ ( S I N ( BETA2-THET2 ) ) 

P3=P2* ( 1+ ( 2 . 8/2 . 4 ) * ( MNN2* *2-1 ) ) 

P03=P3*( ( 1+0.2*M3**2) ) **3.5 

THETT2=THET2/ XF’I 

write (6 .,70) THETT2 , MN3 , M3 

■format (IX, ' THET2 = ' , F6 . 2 , 3X , ' MN3 = ' , F6 . 3, 3X , ' M3= ' , F6 . 3 ) 
read ( 5 , * ) 

X2=XL2*C0S ( THET2+THET1 ) 

Y2=XL2*S IN ( THET2+THET1 ) 
write (6, 61) X2,Y2 

f ormat ( IOX , ' X2= ' , F6 . 3 , 3X 5 ' Y2= ' , F6 . 3 ) 

BETA3=AT AN ( ( Y-Y1-Y2 ) /( X-X1-X2 ) ) -THET1-THET2 
A2=( ( 2/ TAN ( BET A3 ) * ( M3* *2 . * ( S IN ( BET A3 ) ) **2-1 . ) ) ) 

THET3=AT AN ( A2/ ( M3* *2 . * ( 1 . 4+COS ( 2* BET A3 ) ) +2 ) ) 

MNN3=M3*S I N ( BET A3 ) 

MN4=SQRT( (MNN3**2+5)/(7*MNN3**2-l) ) 

M4=MN4/ (SIN( BETA3— THET3) ) 

P4=P3* ( 1+ ( 2 . 8/2 . 4 ) * ( MNN3* *2-1 ) ) 

P04=P4* ( ( 1+0 . 2*M4**2 ) )**3.5 
THETT3=THET3/XPI 


write (6 ,80) THETT3 , MN4 , M4 

format (ix, ' THET3 = ' , F6 . 2 , 3X , ' MN4 = ' , F6 . 3 , 3X , ' M4= ' , F6 . 3 ) 
BETA3-BETA3/XPI 
write ( 6 , * ) ' BETAS = ' , BETAS. 

BETA3=BETA3*XPI 
read ( 5 , * ) 

X3=XL3*CQS ( THET3+THET2+THET 1 ) 

Y3=XL3*S IN ( THET3+THET1+THET2 ) 
write [6 ,62) X3,Y3 

f ormat ( IOX , ' X3= ' , F6 . 3, 3X , ' Y3= ' , F6 . 3) 

BETA4=AT AN ( ( Y-Y1-Y2-Y3 ) / ( X-X1-X2-X3) ) -THET 1-THET2-THET3 
A3=( (2/TAN ( BETA4 ) * ( M4* *2 . * ( SI N ( BETA4 ) )**2-l. ) ) ) 

THET 4=AT AN ( A3/ ( M4**2 . * ( 1 . 4+COS ( 2*BETA4 ) ) +2 ) ) 

MNN4=M4 *S I N ( BET A4 ) 

MN5=SQRT ( (MNN4**2+5)/(7*MNN4**2-l ) ) 

M5=MN5/ (SIN( BETA4— THET4 ) ) 

P5=P4* ( 1+ ( 2 . 8/2 . 4 ) * ( MNN4* *2-1 ) ) 

P05=P5* ( ( 1 +0 . 2*M5 * *2 ) ) * *3 . 5 r 

THETT4=THET4/XPI 
write (6, 90) THETT 4 , MN5 , M5 

format (ix, ' THET 4 =’ ' ,F6.2 ! ,3X, ' MN5 =' ,F6.3,3X, 'M5=' ,F6.3) 
BET A4=BET A4/ XF' I 
write ( 6 , * ) ' BETA4 = ' , BET A4 

BET A4-BET A4*XP I 
read ( 5, * ) 

X4=XL4*C0S ( THET4+THET1+THET2+THET3 ) 

Y4=XL4*S IN ( THET4+THET 1+THET2+THET3 ) 
write(6 ? 63) X4 5 Y4 




BETA5=AT AN ( ( Y-Y1-V2-Y3-Y4 ) / ( X-X1-X2-X3-X4 ) ) 

BETA5=BET A 5-THET 1 — THET2— THET3— THET 4 

A4=( (2/TAN(BETA5}*(M5#*2.*(SIN(BETA5> )**2-l. ) ) ) 

THET5=ATAN(A4/(M5**2.*(1.4+C0S(2*BETA5) ) +2 ) ) 

MNN5=M5#SIN ( BETAS ) 

MN6=SQRT( (MNN5##2+5)/(7*MNN5**2-l) ) 

M6=MNA/(SIN(BETA5-THET5) ) 

P6=P5* ( 1+ ( 2 . 3/2 . 4 ) * ( MNN5**2-i ) ) 

PQ6=P6#< (1+0„2*M6**2) )**3.5 

THETT5=THET5/XP I 

write ( £> , 91 ) THETT5 , MN6 , MS 

format (IX, ' THET5 = ' ,F6.2,3X, 'MN6 = ' ,F6 .3, 3X , ' M6= ' , F6 , 3) 

BETA5=BETA5/XPI 

write(6,*) ' BETAS = ' , BET A5 

BETA5=BETA5*XPI 

T0TPD=P06/P01 

PD=P6/P1 

write (6,*) ' THE TOTAL PRESSURE DRPUT QBLIQi .FOROP IS ' , TOTF'D 

write(6,#) ' THE PRESSURE RISE IS ' ,PD 

read ( 5 , % ) 

stop 

end 



MACH # 4 


program INLETLENGTH 
dimension M, BET, THET 

rea 1 Ml , MN 1 , M2 , MN2 , M3 , MN3 , M4 , MN4 , M5 , MN5 , M6 , MN6 , M7 , MN7 , MNN 
rea 1 MNN 3 , MNN4 , MNN5 , MNN6 , MNN7 
XF'I=3. 141593/180. 

write (6,,#) ' ENTER THE VALUES FOR THETA1 AND Ml' 

read (5,*) THETA, Ml 

write (6,*) 'ENTER THE VALUE OF Y AND L1-L5 ' 

read ( 5 , * ) Y , XL1 , XL2 , XL3 , XL4 , XL5 

Ml=3 . 0 

THETA=8 . 00 

Y=3 . 5 

XL1=7. 031008 
XL2=2. 04934 137 
L3= 1.70362 
L4=l . 10302 
XL3=1 .09605 
XL4=1 .10302 
Pi=i.O 

BET A=THET A+3 . 

THET1=THETA*XPI 
BET1=BETA*XPI 
TT=T AN ( THET 1 ) / 2 . 

TB=0 . 

if (ABS(TT-TB) . It. .00001 ) goto 20 
BET 1=BET1+ ( TT-TB ) *1 . 7 

A=(M1**2*(1.4+C0S(2.*BET1) )+2. )*TAN(BET1) 

XXX=(SIN( BET 1 ) )**2 

TB= ( Ml % ( (SIN(BETl) )*#2)“1. )/A 

goto 10 

write ( 6 , * ) ' BET 1 =',BET1 

MN1=M1*SIN(BET1 ) 

write (6,*) ' MN1 = ' ,MN1 

MN2=SQRT ( (MNl**2+5)/(7*MNl**2-l) ) 

write ( 6 , * ) ' MN2 =',MN2 

write ( 6 , * ) BET 1 , THET 1 

M2=MN2/ (SIN ( BET 1 -THET 1 ) ) 

P2=Pl*(l+(2.3/2.4)*(MNl**2-l) ) 

P01=P1* ( ( 1+0.2*M1**2) ) **3.5 
P02=P2* ( ( 1+0 . 2*M2**2 ) )**3.5 
write ( 6 , * ) M2 
continue 
BETA=BET 1 / XPI 

write ( 6 , * ) ' BETA = ' ,BETA 

BETA=BETA*XPI 

X=Y/TAN ( BETA ) 

X1=XL1 *C0S ( THET 1 ) 

Y1=XL1*SIN( THET 1 ) 

BETA2=ATAN( ( Y-Yl ) / ( X-Xl ) ) -THET 1 

BETAI-BETA2/XPI 

write (6, 50) X,Y 

format ( ioX ,' X = ' , F6 . 3 , 3X Y = ',F6.3) 
wr i te ( 6 , 60 ) BETA2,X1,Y1 

format ( ioX , ' BET A2= ' ,F6.2,3X , ' Xl=' ,F6.3,3X , ' Yl=' „F6.3) 

BET A2=BET A2 * XPI 
read (5,*) 

Al=( ( 2/TAN ( BETA2 ) * ( M2* *2 . * ( SIN ( BETA2 ) )**2-l. ) ) ) 

THET2=AT AN ( Al/ ( M2* *2 . * ( 1 . 4+COS ( 2* BET A2 ) ) +2 ) ) ' 


write (6,*)' MNN2 = , MNN2 

MN3=SGRT( (MNN2**2+5)/(7*MNN2**2-l) ) 

M3=MN3/ ( S I N ( BETA2-THET2 ) ) 

P3*P2*(l+(2.8/2.4)#(MNN2**2-l ) ) 

PQ3=P3*( ( 1+0.2*M3**2) )**3.5 

THETT2=THET2/XPI 

write (6, 70) THETT2 , MN3 , M3 

format ( IX , ' THET2 = ' , F6 . 2 , 3X , ' MN3 = ' ,F6 . 3 , 3X , ' M3= ' , F6 . 3 ) 
read ( 5 , * ) 

X2=XL2*C0S(THET2+THET1 ) 

Y2=XL2*SIN(THET2+THET1 ) 
write ( 6 , 61 ) X2,Y2 

format ( iOX , ' X2= ' , F6 . 3, 3X , ' Y2= ' , F6 . 3 ) 

BETA3=ATAN ( ( Y-Y1-Y2 ) / ( X-X1-X2 ) ) -THET1-THET2 
A2= ( ( 2/TAN (BETA3)*(M3**2.*< SIN (BETA3) )**2-i. ) ) ) 
THET3=ATAN < A2/ ( M3**2 . * ( i . 4+COS ( 2*BETA3 ) ) +2 ) ) 

MNN3=M3*S I N ( BET A3 ) 

MN4=SQRT( (MNN3**2+5)/(7*MNN3**2-l) ) 

M4=MN4/ ( SI N ( BET A3-THET3 ) ) 

P4=P3* ( 1+ ( 2 . 8/2 . 4 ) * ( MNN3* *2-1 ) ) 

P04=P4* ( (i+0„2*M4**2) )**3.5 
THETT3=THET3/XPI 

write (6, 30) THETT3 , h1N4 , M4 

format ( IX , ' THET3 = ' , F6 . 2 , 3X , ' MN4 = ' ,F6. 3,3X , ' M4= ' ,F6 .3) 

BETA3=BETA3/XPI 

write(&,*) ' BETAS = ',BETA3 

BET A3=BET A3* X P I 

read ( 5 , * ) 

T0TPD=P04/P01 

PD=P4/P1 

wri te ( 6 , * ) ' THE TOTAL PRESSURE DROP IS , TOTPD 

write ( 6 , * ) ' THE PRESSURE RISE IS ' ,PD 

read ( 5 , * ) 

stop 

end 



MACH # 3 AND MACH # 2 


program INLETLENGTH 
dimension M, BET, THET 

rea 1 Mi , MN1 , M2 , MN2 , M3 , MN3 , M4 , MN4 , M5 , MN5 , M6 , MN6 , M7 , MN7 , MNN2 
rea 1 MNN3 , MNN4 , MNN5 MNN6 , MNN7 
XPI=3. 141593/180. 

write (6,*) ' ENTER THE VALUES FDR THETA1 AND Ml' 

read (5,*) THETA, Mi 

write ( 6 ,%) 'ENTER THE VALUE OF Y AND L1-L5 ' 

read (5 , * ) Y , XL1 , XL2 , XL3 , XL4 , XL5 

Ml =3.0 

THETA=3 . 00 

Y=3 . 5 

FOR MACH # 3 
XL 1=6 . 0 


FOR MACH 2 
XLi=3 . 0 

XL 1=9. 13034937 
XL2=2. 04934137 
L3= 1.70362 
L4=l . 10302 
XL3= 1.09605 
XL4=1. 10302 
Pl=l .0 

BET A=THET A+3 . 

THET1=THETA*XPI 
BET1=BETA*XPI 
TT =TAN ( THET 1 ) /2 . 

TB=0 . 

if (ABS(TT-TB) . It. .00001) goto 20 
BET 1=BET 1+ ( TT-TB ) *1 . 7 

A=(M1**2#(1.4+C0S(2.*BET1) )+2. )*TAN(BET1) 

XXX= ( S IN ( BET 1 ) ) **2 

TB=(M1**2*( ( SI N ( BET 1 ) )**2)-l. ) /A 

goto 10 

write(6,*) ' BET 1 =' ,BETi 
MN1=M1*SIN(BET1 ) 
write(6,*) ' MN1 = ',MN1 
MN2=SQRT( (MNl**2+5)/(7*MNl**2-l) ) 
write (6,*) ' MN2 =',MN2 
write ( 6 , * ) BET 1 , THET 1 
M2=MN2/(SIN(BET1-THET1) ) 

P2=P1 * ( 1+ ( 2 . 8/2 . 4 ) * ( MN1*#2-1 ) ) 

P01=P1*( (1+0.2#M1*#2) ) **3 . 5 
P02=P2* ( < 1+0.2*M2#*2) ) **3. 5 
write ( 6 , ) M2 
continue 
BETA=BET 1 / XP I 

write ( 6 , tfc ) ' BET A = '.BETA 

BET A=BETA*XP I 
X=Y/TAN ( BETA ) 

X1=XL1*C0S(THET1 ) 

Y1=XL1*SIN(THET1 ) 

BETA2=ATAN ( ( Y-Yl ) / ( X-Xl ) ) -THET 1 
BET A2=BET A2 / X P I 
write (6, 50) X,Y 



write (6, 60) BETA2,X1,Y1 

format ( iox , ' BET A2= ' , F6 . 2 , 3X , ' X 1= ' , F6 . 3 , 3X , ' Yl= ' , F6 . 3 ) 

BETA2=BETA2*XPI 

read ( 5, * ) 

Al=( (2/TAN(BETA2)*(M2**2.*(SIN(BETA2) )#*2-l. ) ) ) 
THET2=ATAN(Ai/(M2**2.*(1.4+C0S<2#BETA2) )+2) ) 

MNN2=M2*S I N ( BETA2 ) 

write (6, % ) ' MNN2 = ' , MNN2 

MN3=SQRT ( (MNN2**2+5)/(7*MNN2**2-i) ) 

M3=MN3/ ( SIN ( BET A2-THET2 ) ) 

P3=P2* ( 1+ ( 2 . 8/2 . 4 ) * ( MNN2**2-1 ) ) 

P03=P3*( (1+0.2*M3##2) )**3.5 

THETT2=THET2/XPI 

write (6, 70) THETT2 , MN3 , M3 

format { ix , ' THET2 = ' ,F6.2,3X , ' MN3 =' ,F6.3,3X , ' M3=' ,F6.3) 
read ( 5 , * ) 

T0TPD=P03/P01 

PD-P3/P1 

write ( 6 , * ) ' THE TOTAL PRESSURE DROP IS ',T0TPD 

write ( 6 , * ) ' THE PRESSURE RISE IS ' ,PD ! 

read ( 5 , * ) 

stop 

end 





0 
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il 


10 


real ‘LAMBD , MHU , N , MACH I , LAEF , LEF ^ ~ 

RH0IN=0 . 0000321 14 
RHOAT=0 . 002367 
' VELIN=9S13.2 
RADN=0.5 
LAMBD=69 . 53 
ALPHA=0 . 0 
AN0FA=0 . 0 
MHU=0 . 0000093 
N=0 . 00 
MACH I =10.0 
BETA=0 . 0 
RADLE=0 . 02083 
EMISS=0.S9 
STEFB= . 000000001714 
BAMMA=1 . 34 
TN0T=753S . 2 

HOLAM= 0 . 1055* ( RHO I N/RHOAT ) * * . 5* ( VEL I N/ 1 0* *4 ) * * 1 . 16 
H0TUR=0.437* (RHOIN/RHOAT) ** .78* (VELIN/10**4 ) **1 . 54 
HNOSE-HOLAM/RADN** . 5 
LAEF=S I N ( LAMBD ) *COS ( ALPHA } 

LEF=AS I N ( LAEF ) 

PF1= ( 1 . 33#MACH I * *2+2 . 5 ) * ( ( COS ( LEF ) ) * *2+0 . 0019 ) 

PF2=1 . 33* ( MACH I * ( COS ( LEF ) ) )##2+l 
PF=PF2/PF1 


HL1=H0LAM* ( PF/RADLE ) ** . 5 

HL2=0 . 72 * ( COS ( LEF ) ) * * 1 . 1 5+ . 04 *S I N ( LEF ) 

HLSTAG=HL.1*HL2 
dss 10 J0=1 , 106 
XX=J0/2 .0 
XDIS=J0*1 .0 

RN=RHOIN*VELIN*XDIS/MHU 
if ( RN . 1 t . 3 . 5E5 ) then 

HLAM1=H0LAM* (PF*(1+2*N)/(3.*XDIS))**.5 
HLAM2=1 .75*SIN( BETA ) - . 86* ( ( SIN ( BETA ) )**2> 

HL AM3=0 . 036/ ( ( MACH I *S I N ( BETA > ) * *2+1 . 0 ) 

HL AM=HLAM1 * ( HLAM2+HLAM3 ) 

UE=VELIN*SIN ( LEF ) 

TAW=TNOT— ( 1— . 15 ) *UE* *2/ 12012 . 0 
goto 20 
endif 

HTURB1=H0TUR# ( XDI S/10 . 0 ) * * ( -2 . 0 ) * ( ( 1 . 0+1 . 25*N > / 
HTURB2=PF * * O . 8 * ( 2 . 9 * ( S I N ( 1 . 5 * BET A ) ) - 1 . 6 # ( ( S I N ( B 
HTURB3=0 . 02 / ( ( MACH I * 3 1 N ( BET A ) ) * * 2+ 1 . O ) 
HTURB=HTURB1 * ( HTURB2+HTURB3 ) 

TAW=TNOT— 0. 1*VELIN**2* ( (SIN (LEF) ) **2 ) / 12012 . 0 

HWU=HTURB 

goto 80 

HWU=HLAM 

do 90 J=1 , 1E4 , 10 
TT=J*1 . 0 

A=TT * *4+HWU*TT / ( EMI SS*STEFB > 
B=HWU*TAW/(EMISS*STEFB) 

AB=A— B 

i f ( AB . g t . 1 . OE -6 ) t hen 
goto 110 
endif 

if ( AB . 1 t . — 1 . OE— 6 ) then 
TT =TT +10 . O 
goto 111 
endif 
goto 40 
continue 

write (6,100) HWU , TT*20 , XDI S/2 
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LIST OF ENERGY STATE FORTRAN PROGRAMS 

Bob Stonebraker 


PROGRAM HECONST - Computes lines of constant He (ft.)* 

PROGRAM FSCONST - Computes contours of constant Fs (ft. /lb.) 

for phase 1 of the mission. 

PROGRAM FSCNST2 - Computes contours of constant Fs (ft. /lb.) 

for Phase 2 of the mission. 

PROGRAM PHASE1 - Computes fuel consumption, elapsed time and 

range for Phase 1. 

PROGRAM PHASE2 - Computes fuel consumption, elapsed time and 

range for Phase 2 . 

SUBROUTINE ATMOS FR (ALT,T,RO,P) - Computes standard atmospheric 

variables as a function of altitude. 

SUBROUTINE TFRMJT (M, ALT, FN, SFC) - Computes net thrust and sfc 

from given data as a function of Mach number 
and altitude. (based on economic sfc) 

SUBROUTINE SCRMJT (M, ALT, FN, SFC) - Computes net thrust and sfc 

from given data as a function of Mach number 
and altitude. (based on equivalent fuel/air 
ratio of 0.8) 

SUBROUTINE DRAG (M, ALT, D) - Computes drag from given data as a 

function of Mach number and altitude. 

SUBROUTINE SPLINE (X, F, XF, FN,N) - Interpolates between points 

(X, F) for point (XF, FN) . 
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PROGRAM HECONST 
ROBERT STONEBRAKER 

COMPUTES VALUES OF CONSTANT He (ft.) 
IMPLICIT REAL ( A-H , O-Z ) 

REAL M 

OPEN ( 5 , FILE= ' HECONST . DAT ' ) 

R=1715 .621 
G=32 . 1578 
HE-2 . E4 
DO 200 1=1,10 
HG=0 . 

100 CALL ATMOSFR(HG,T,RO,P) 

M=SQRT ( (HE-HG) *2 . *G/ ( 1 . 4*R*T) ) 

WRITE ( 5 , * ) M,HG,HE 
WRITE ( * , * ) M,HG,HE 
HG=HG+4000 . 

IF(HG.GT. 110000.) GOTO 200 
IF(HG.LT.HE) GOTO 100 
200 HE=HE+2 . 0E4 
CONTINUE 
STOP 
END 


oonoo 


PROGRAM FSCONST 
ROBERT STONEBRAKER 

COMPUTES VALUES OF CONSTANT Fs (ft/lb) OVER PHASE 1 
WITH NET THRUST, SFC, AND DRAG AS FUNCTIONS OF 
ALTITUDE AND MACH NO. 

IMPLICIT REAL(A-H,0-Z) 

REAL M 

OPEN ( 5 , FILE= ' FSCONST . DAT ' ) 

WRITE (*,*) 'ENTER INITIAL VEHICLE WIEGHT IN LBS.' 

READ ( * , * ) W 
M= . 8 

R=1715 .621 
C MACH NO. LOOP 

DO 100 1=1,27 
ALT=40000. 

C ALTITUDE LOOP 

DO 200 J=l, 30 

CALL ATMOS FR ( ALT , T , RO , P ) 

A=3600*SQRT ( 1 . 4*R*T) 

C 

CALL TFRMJT (ALT,M, FN , SFC) 

IF(FN.EQ. 0. ) GOTO 200 
IF (SFC. EQ. 0. ) GOTO 200 
C 

CALL DRAG (ALT, M,D) 

FS=M*A* ( FN-D) / ( FN *SFC*W) 

IF (FS . LT. 0 . ) THEN 
FS=0 . 

ENDIF 

WRITE (* , 10) M , ALT , FN , D , SFC , FS 
10 FORMAT (F5 . 2 , 3X, F8 . 1 , 3X, F10 . 1 , 3X, F10 . 1, 3X, F6 . 3 , 3X, F7 . 1) 
WRITE (5,20) M , ALT , FS 
20 FORMAT (F5 . 2 , 3X, F8 . 0, 3X, F7 . 1) 

200 ALT=ALT+2000 . 

W=W-100. 

100 M=M+. 200000 
STOP 
END 
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PROGRAM FSCNST2 
ROBERT STONEBRAKER 

COMPUTES VALUES OF CONSTANT Fs (ft/lb) OVER PHASE 2 
WITH NET THRUST, SFC, AND DRAG AS FUNCTIONS OF 
ALTITUDE AND MACH NO. 

IMPLICIT REAL ( A-H , O-Z ) 

REAL M 

OPEN ( 5 , FILE= ' FSCNST2 . DAT • ) 

WRITE (*,*) 'ENTER INITIAL VEHICLE WIEGHT IN LBS.' 

READ ( * , * ) W 
M=6 . 

R=1715 .621 
MACH NO. LOOP 
DO 100 1=1,41 
ALT=70000. 

ALTITUDE LOOP 
DO 200 J=1 , 16 
SPEED OF SOUND (ft/sec) 

CALL ATMOSFR ( ALT , T , RO , P) 

A=SQRT (1.4 *R*T) 

Q=. 5*RO*l. 4*R*T*M**2 
IF(Q. GE . 1800) THEN 
FS=0 

GOTO 300 
ENDIF 

CALL SCRMJT (ALT,M, FN, SFC) 

FN=4*FN 
SFC=SFC/ 3600 
IF (FN. EQ. 0 . ) GOTO 200 

CALL DRAG (ALT, M,D) 

FS=M*A* (FN-D)/ (FN*SFC*W) 

WRITE (*,10) M, ALT, FN , D, SFC, FS 

FORMAT (F4 . 1, 3X, F7 . 0 , 3X, F10 . 1 , 3X, F10 . 1 , 3X, F9 . 6 , 3X, E12 .5) 
WRITE (5,20) M, ALT, FS 
FORMAT (F4 . 1 , 3X, F7 . 0, 3X, F12 . 3) 

ALT=ALT+ 2000. 

W=W-147. 

M=M+ . 10 
STOP 
END 
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PROGRAM PHASE1 
ROBERT STONEBRAKER 

COMPUTES THE FUEL CONSUMED ALONG PREDETERMINED FLIGHT PATH 
FROM (M=. 8 AT 40,000' TO M=6 AT 76,000') 

IMPLICIT REAL(A-H,0-Z) 

REAL M 

OPEN ( 5 , FILE= ' PHASE 1 . DAT ' ) 

OPEN (6 , FILE= ' PHASE 11 . DAT ' ) 

WRITE(*,*) 'ENTER INITIAL VEHICLE WEIGHT' 

READ ( ★ 9 ★ ) W 

WRITE (*,*) 'ENTER MACH NO. STEP SIZE* 

READ ( * , * ) DM 
R=1715 .621 
G=32 . 1578 
M= . 8 

HG=40000. 

HGO=40000. 

A=968 . 1 

C COMPUTE INITIAL HE AND F=1/FS AT M=.8 AND ALT (HG) =40000 
HE1=49326 . 186 
CALL TFRMJT (HG,M, FN, SFC) 

SFC=SFC/3600. 

CALL DRAG (HG,M, D) 

P1=W/(M*A*(FN-D) ) 

F1=FN*SFC*P1 

C 

N=5 . 2/DM+l 
DO 100 1=1, N 
M=M+DM 

CALL FLTPTH(M,HG) 

CALL ATMOSFR (HG , T , RO , P) 

CALL TFRMJT (HG,M,FN, SFC) 

CALL DRAG(HG,M, D) 

SFC=SFC/3600 . 

A=SQRT (1.4 *R*T) 

HE2=HG+.5*(1.4*R*T*M**2)/G 

TR=FN-D 

IF(TR.EQ.O) GOTO 100 
P2=W/ (M*A* (FN-D) ) 

F2=FN*SFC*P2 

DTIME= (HE2-HE1) * .5* (P2+P1) 

DFUEL= (HE2-HE1) * . 5* (F2+F1) 

W=W-DFUEL 
TIME=TIME+DTIME 
ACC=DM/ (M*DTIME) 

RNG=RNG+SQRT ( (M*A*DTIME) **2- (HG-HGO) **2 ) /5280 . 

HGO=HG 

FUEL=FUEL+DFUEL 

ELT=TIME/60 

WRITE ( * , 10) ELT , RNG , M , HG , FN , D , W 
WRITE (5, 10) ELT , RNG , M , HG , FN , D , W 
10FORMAT (F7 . 3 , 2X, F8 . 1 , 2X, F6 . 2 , 2X, F8 . 0 , 2X, F8 . 0 , 2X, F8 . 0 , 2X, F8 . 0) 
WRITE (6,20) RNG, M,HG, ACC, FUEL 
20 FORMAT (F8 ,1,2X,F6.2,2X,F8.0,2X,F7.4,2X,F7.0) 




HE1=HE2 
P1=P2 
100 F1=F2 

TIME=TIME/60 
WRITE ( * , * ) 'M 
WRITE ( * , * ) 'T 


“'.H,' ALT 
>* ' ,TIME, 'min. 




=b • » ffr l 

*W fuel =', FUEL, 'lbs. • 


\s 1 

\ 



o o o o o o o 


PROGRAM PHASE 2 
ROBERT STONEBRAKER 

COMPUTES FUEL CONSUMED ALONG PREDETERMINED FLIGHT PATH 
FROM (M=6 AT 76,000* TO M=10 AT 100,000*) 

IMPLICIT REAL (A-H,0-Z) 

REAL M 

OPEN ( 5 , FILE= ' PHASE2 . DAT » ) 

OPEN ( 6 , FILE= ' PHASE22 . DAT • ) 

WRITE (*,*) ‘ENTER INITIAL VEHICLE WEIGHT* 

READ ( * , * ) W 

WRITE(*,*) 'ENTER MACH NO. STEP SIZE* 

READ ( * , * ) DM 

WRITE (*,*) 'ENTER NO. OF SCRAMJETS ' 

READ ( * , * ) K 
R=1715 . 621 
G=32 .1578 
M=6 

HG=75835. 

HGO=75835. 

A=968 . 1 

COMPUTE INITIAL HE AND F=1/FS AT M=6 
HE1=600433. 

CALL SCRMJT(HG,M, FN,SFC) 

FN=K*FN 
SFC=SFC/3600. 

CALL DRAG(HG,M, D) 

P1=W/(M*A*(FN-D) ) 

F1=FN*SFC*P1 

N=4 ./DM+30 
DO 100 1=1, N 
M=M+DM 

FLIGHT PATH ALONG Q=1800 
IF (M.GE. 10 . ) THEN 
M=10. 

HG=HG+300. 

IF(HG.GE. 100000. ) GOTO 500 
GOTO 50 
ENDIF 

HG=-7.35005*M**4+262.753*M**3-3749.61*M**2+29936.5*M-16027. 
C 

50 CALL ATMOSFR (HG , T , RO , P) 

CALL SCRMJT (HG,M, FN , SFC) 

FN=K*FN 

SFC=SFC/3600. 

CALL DRAG(HG,M,D) 

A=SQRT (1.4 *R*T) 

HE2=HG+ . 5* ( 1 . 4 *R*T*M**2 ) /G 
TR=FN-D 

IF(TR.EQ.O) GOTO 100 
P2=W/(M*A*(FN-D) ) 

F2=FN*SFC*P2 

DTIME= (HE2-HE1) * . 5* (P2+P1) 





DFUEL* (HE2-HE1) * . 5* (F2+F1) 

W=W-DFUEL 
TIME=TIME+DTIME 
ACC=DM/ (M*DTIME) 

RNG=RNG+SQRT( (M*A*DTIME) **2-(HG-HG0) **2)/5280. 

HGO=HG 

ELT=TIME/60 

FUEL=FUEL+DFUEL 

WRITE (*, 10) ELT, RNG,M, HG, FN, D, W 
WRITE (5, 10) ELT , RNG , M , HG , FN , D , W 
10FORMAT (F7 . 3 , 2X, F8 . 1 , 2X, F6 . 2 , 2X, F8 . 0 , 2X, F8 .0, 2X, F8 . 0 , 2X, F8 . 0) 
WRITE (6, 20) RNG,M,HG, ACC, FUEL 
20 FORMAT (F7 .3,2X,F6.2, 2X, F8 . 0, 2X, F7 . 4 , 2X, F7 . 0) 

HE1=HE2 
P1=P2 
100 F1=F2 

500 WRITE (*,*) 'FUEL =' , FUEL 
STOP 
END 
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SUBROUTINE ATMOS FR (HG , T , RO , P) 

ROBERT STONEBRAKER 

COMPUTES TEMPERATURE, DENSITY, AND PRESSURE GIVEN A 
GEOMETRIC ALTITUDE HG. 

IMPLICIT REAL(A-H,0-Z) 

RE=2 . 086466E7 
H=HG*RE/ (RE+HG) 

IF (H.LT. 36000. )THEN 
H0=0. 

T0=518 . 69 / 

R00=2 . 3769E-3 
P0=2116. 2 
A=-3 . 56627E-3 
T=TEMP (TO, A, H, HO) 

RO=ROS LOP ( ROO , T , TO , A) 

P=PSLOP ( PO , T , TO , A) 

GOTO 500 
ENDIF 

IF (H.LT. 82000.) THEN 
H0=36000 
T=389 . 99 
R00=7 . 0858E-4 
P0=474 .7098 
RO=ROLNR ( ROO , T , H , HO ) 

P=PLNR(P0,T,H,H0) 

GOTO 500 
ENDIF 

IF (H.LT. 154000.) THEN 
H0=82000 
T0=389 . 99 
R00=7 . 7664E-5 
P0=52 . 03 
A=l. 64355E-3 
T=TEMP (TO, A, H, HO) 

RO=ROSLOP (ROO , T , TO , A) 

P=PSLOP ( PO , T , TO , A) 

GOTO 500 
ENDIF 

500 RETURN 
END 
C 

FUNCTION TEMP (TO, A, H, HO) 

TEMP=T0+A* (H-HO) 

RETURN 

END 

FUNCTION ROSLOP (ROO , T , TO , A) 

ROS LOP=RO 0 * ( (T/TO) ** (- (32 . 1578/ (A*1715 . 621) +1 . ) ) ) 

RETURN 

END 

FUNCTION PSLOP(PO,T,TO,A) 

PSLOP=PO* ( (T/TO) ** (-32 . 1578/ (A*1715. 621) ) ) 
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RETURN 

END 

FUNCTION ROLNR(ROO,T,H,HO) 

ROLNR=ROO*EXP (“32 . 1578/(1715 . 621*T) * (H-HO) ) 

RETURN 

END 

FUNCTION PLNR ( PO , T , H , HO ) 

PLNR=PO*EXP (-32. 1578/ ( 1715. 621*T) * (H-HO) ) 

RETURN 

END 



SUBROUTINE TFRMJT(HG,M,FN,SFC) 

C ROBERT STONEBRAKER 

C COMPUTES THE NET THRUST AND SFC OF THE FULL SCALE 
C TURBOFAN-RAMJET AS A FUNCTION OF ALTITUDE AND MACH NO. 

C BASED ON CURVE FIT POLYNOMIALS AND NATURAL CUBIC SPLINES. 

C 

IMPLICIT REAL (A-H,0-Z) 

REAL M 

DIMENSION HGA(15) ,FNA(15) ,SFCA(15) 

C 

IF(M.LT. 0.8) THEN 
FN=0 
SFC=0 
GOTO 500 
ENDIF 

IF (M. GT. 6 . 0) THEN 
FN=0 
SFC=0 
GOTO 500 
ENDIF 
C 

K=7 

C *** 0.8 TO 3.5 
HGA(l) =40000 

FNA (1) =-2879 . 52*M**5+28051 . 4*M**4-101019 . *M**3+172092 . *M**2 

* -112116. *M+37378. 7 

SFCA(l) =-. 0430727*M**6+. 570346*M**5-3 . 00762*M**4+8 . 02043*M**3 

* -11. 2762*M**2+7.86021*M-1. 40087 
C 

C *** 2.0 TO 4.5 
HGA(2) =50000 

FNA(2) =2569 . 3 3 *M* *3 -152 06 . 7*M**2+53954 . 5*M-37169 . 1 
IF (M. LE. 3 . 5) THEN 

SFCA (2 ) = . 0640001*M**3-. 492*M**2+1 . 3 16 *M-. 424001 
GOTO 60 
ENDIF ' 

SFCA (2 ) =. 034*M**2-. 237*M+1. 312 
C 

C *** 2.5 TO 6.0 
60 HGA(3) =60000 

FNA (3)=-905. 152 *M** 4+14 058 . 4*M**3-72866 . 3*M**2+174728 . *M-139344 . 
IF (M. LT. 3.5) THEN 

SFCA(3 ) =. 0880001*M**2-. 422*M+1 .299 

GOTO 70 

ENDIF 

SFCA ( 3 ) = -.0020003*M**3+.0522861*M**2-.28693*M+1. 34922 
C 

C *** 3.0 TO 6.0 
70 HGA (4 ) =70000 

FNA ( 4 ) =-1287 . 11*M**4+22021 . 3*M**3-134968 . *M** 2+3 7 18 2 4 . *M-369131 . 
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SFCA (4) =. 0103333*M**3-. 13 092 9 *M* *2+ . 6 13 02 4 *M- . 1077 

*** 4.0 TO 6.0 
HGA(5) =80000 

FNA(5) =-2007 . 67*M**3+30631 . 5*M**2-137628 . *M+212654 . 
SFCA(5) =. 00866667*M**3-. 105 *M**2+. 48 6334 *M+. 0959986 

*** 4.0 TO 6.0 
HGA(6) =90000 

FNA(6) =-1233 . 33*M**3+18722 . *M* *2 -8 3 600. 7*M+128444 . 
SFCA(6) =. 0076661*M**3-. 0854991*M**2+. 37282 9*M+. 311008 

*** 4.0 TO 6.0 
HGA (7 >=100000 

FNA(7) =-683 . 667*M**3+10277 . 5*M**2-44661 . 8*M+66948 . 
SFCA(7)=.024*M**3-.320001*M**2+1.476*M-1. 37301 

CALL SPLINE (HGA, FNA,HG,FN,K) 

IF (FN. LT. 0 . ) THEN 
FN=0 
ENDIF 

CALL SPLINE (HGA, SFCA, HG, SFC, K) 

IF (SFC. LT. 0 . ) THEN 
SFC=0 
ENDIF 
500 RETURN 
END 
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SUBROUTINE SCRMJT(HG,M,FN,SFC) 

ROBERT STONEBRAKER 

COMPUTES FG(lbsf), FN(lbsf), ISP (sec) , SFC(lbsm/lbsf/hr) 

AS FUNCTION OF ALT. AND MACH NO. 

CURVE FIT EQUATIONS BASED ON FUEL/AIR RATIO OF 0.8 
IMPLICIT REAL (A-H,0-Z) 

IMPLICIT REAL (M,I) 

DIMENSION MA(15) ,F(15) 

DATA MA(1) ,MA (2) ,MA(3) /5 .0,7.5,10.0/ 

N=3 

*** pg *** 

F (1) =-4 . 92613E-10*HG**3+1. 45585E-4*HG**2-14 . 8904*HG+534062 
F (2 ) =— 4 . 91975E-10*HG**3+1 . 71206E-4*HG**2— 20 . 4451*HG+845486 
F(3)=-4.99582E-10*HG**3+1.95693E-4*HG**2-26.1568*HG+1201470 
CALL SPLINE (MA, F,M, FG,N) 

*** FN *** 

F (1) =-l. 61294E-10*HG**3+4 . 7169 2E-5*HG* *2-4 . 77752*HG+169622 
F(2)=-1.04412E-10*HG**3+3.57416E-5*HG**2-4.20399*HG+171214 
F (3) =-4 . 13576E-11*HG**3+1 . 7587E-5*HG**2-2 . 52803 *HG+123291 
CALL SPLINE (MA, F,M, FN,N) 

*** ISP *** 

F(1)=-5.79888E-12*HG**3+1.19854E-6*HG**2-8.53466E-2*HG+5398.13 

F(2)=-9.88777E-12*HG**3+2.55155E-6*HG**2-0.22344*HG+9142.38 
F(3) =-7 . 61882E-12*HG**3+1 . 98524E-6*HG**2— 0 . 174684*HG+7093 . 69 
CALL SPLINE (MA, F,M, ISP,N) 

C 

SFC=3600 ./ISP 

RETURN 

END 
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SUBROUTINE DRAG(HG,M, D) 

C ROBERT STONEBRAKER 

C COMPUTES DRAG OF 515H VEHICLE DESIGN AS A FUNCTION OF 

ALTITUDE 

C AND MACH NO. BASED ON CURVE FIT POLYNOMIALS OF GIVEN DATA. 
IMPLICIT REAL(A-H,0-Z) 

REAL M 

DIMENSION HGAA ( 15 ) , DA ( 15) 

C 

K=13 

HGAA (1) =40000 

DA ( 1) =981 . 718*M**2+10219 . 2*M-5885 .21 
HGAA(2) =45000 

DA(2) =777 . 411*M**2+7952 . 26*M-4088 . 61 
HGAA (3) =50000 

DA (3) =617 . 958*M**2+6144 . 13*M-2528 . 24 
HGAA (4) =55000 

DA(4) =494 . 109*M**2+4690 . 13*M-1114 . 06 
HGAA (5) =60000 

DA (5) =398 . 678*M**2+3505 . 97*M+235 . 57 
HGAA(6) =65000 

DA (6) =326. 16*M**2+2523 . 4*M+1598 . 03 
HGAA (7) =70000 

DA (7) =272 . 36*M**2+1685 . 99*M+3051 .45 
HGAA (8) =75000 

DA (8) =234 . 18*M* *2+945 . 452*M+4678 . 9 
HGAA (9) =80000 

DA (9) =209 . 415*M**2+259 . 422*M+6572 
HGAA (10) =85000 

DA(10) =196 . 64*M**2-409 . 555*M+8834 . 34 
C **NOTE**LAST THREE EQUATIONS ARE NOT VALID FOR M < 4 
IF (M. LT. 4.0) THEN 
K=10 

GOTO 500 
ENDIF 

HGAA(ll) =90000 

DA (11) =-59 . 1967*M**3+1483 . 62 *M**2 -10103 . 1*M+31674 . 9 
HGAA (12) =95000 

DA (12) =-69 . 7467*M**3+1732 . 75*M**2-12589 . 8*M+39067 . 1 
HGAA (13) =100000 

DA (13) =-83 . 6275*M**3+2064 . 89*M**2-15665. 8*M+48294 . 5 
C 

500 CALL SPLINE (HGAA, DA, HG,D,K) 

D=1 . 1*D 

RETURN 

END 
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SUBROUTINE SPLINE (X,F,XN,FX,N) 

ROBERT STONEBRAKER 4/10/90 

COMPUTES VALUES OF AN UNKNOWN FUNCTION FROM UNEVENLY SPACED 
DATA USING A NATURAL CUBIC SPLINE. 

IMPLICIT REAL (A-H,0-Z) 

DIMENSION X(15) ,F(15) ,A(15) ,B(15) ,C(15) ,G(15) ,R(15) 

DETERMINES WHICH INTERVAL OF THE DOMAIN. HOLDS X 
DO 100 1=2, N 
IF(XN.LE.X(I) ) THEN 
J=I-1 
GOTO 200 
ENDIF 

100 CONTINUE 

DEFINES MATRIX COEFFICIENTS 
200 DO 300 1=1, N 

DELX=X(I+1) — X ( I ) 

A(I)=(X(I)-X(I-1) )/DELX 
B(I)=2*(X(I+1) -X(I-l) )/DELX 
C(I)=1 

R(I)=6*(F(I+1)-F(I))/(DELX**2) 

*-6* (F(I) -F (1-1) ) / (DELX* (X (I) -X(I-l) ) ) 

300 CONTINUE 
A(1)=0. 

C(N)=0. 

THOMAS ALGORITHM FOR G(i) 

DO 400 1=2, N 
FACTOR=A(I)/B(I-l) 

B(I)=B(I) -FACTOR*C (1-1) 

400 R(I) =R(I) -FACTOR*R(I-l) 

G(N)=R(N)/B(N) 

DO 500 1=2, N 
NI=N-I+1 

500 G(NI)=(R(NI) -C(NI) *G(NI+1) )/B(NI) 

G(1)=0. 

G (N) =0 . 

DELX=X(J+1)-X(J) 

FX=G ( J) * ( ( (X ( J+l) -XN) **3) /DELX-DELX* (X (J+l) -XN) )/6 
*+G (J+l) * ( ( (XN-X (J) ) **3) /DELX-DELX* (XN-X (J) ) )/6 
*+F ( J) * (X (J+l) -XN) /DELX+F (J+l) * (XN-X ( J) ) /DELX 

RETURN 

END 
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